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Rev. B 58, 7619 (1998)]. By allowing anisotropic interactions in spin space, an 
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universal properties of the models. At the critical point, the massless degrees 
of freedom are described in terms of an effective S=l/2 Heisenberg spin chain 
and two critical Ising models. The Toulouse limit solution enables us to discuss 
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der. Finally, it is shown that the critical point becomes unstable upon switching 
on some weak backscattering perturbations in the frustrated three-leg ladder. 
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I. INTRODUCTION 



One-dimensional quantum spin systems exhibit fascinating physical properties. Most of the 
striking features obsverved in these systems are pure quantum effects and give rise to many 
exotic ground-states or unconventional excitation spectra that are not accessible by traditional 
approaches like spin-wave or perturbation theory. A famous example is the prediction made 
by Haldane M that half-integer antiferromagnetic spin chains are gapless whereas those with 
integer spin have a finite energy gap in its spectrum. 

The different one-dimensional spin gapped phases are usually characterized by the nature 
of the broken discrete symmetries (lattice translational symmetry for instance) in the ground- 
state, the quantum numbers associated with the massive excitations, and the structure of the 
dynamical magnetic susceptibility. Many different spin gapped phases have been identified. 
The spin-f Heisenberg chain is certainly the most famous example of a spin liquid with a gap. 
ft is characterized by a hidden topological order Q related to the breakdown of a hidden Z2 x 
Z2 symmetry ||. The spin excitations are optical magnons carrying spin S=l and display a 
single peak near q — n in the dynamical magnetic susceptibility. The low-lying excitations of 
a massive phase can also have fractional quantum numbers as in the J1-J2 S=l/2 Heisenberg 
chain for a sufficiently large value of J2: The ground-state is spontaneously dimcrized with 
a two-fold degeneracy and the fundamental excitations are massive S=l/2 solitons (spinons) 
separating the two ground-states ||,[3| . Another scenario can also be realized in a two-leg spin 
ladder coupled by a four spin exchange interaction where the elementary excitation is neither 
an optical magnon nor a massive spinon but rather a pair of propagating triplet or singlet 
solitons connecting the two spontaneously dimerized ground-states ^|J7) . This special structure 
of the excitation will reveal itself in the dynamical magnetic susceptibility which displays a 
two-particle threshold instead of a sharp single magnon peak near q = 7r as in the S=l spin 
chain. This incoherent background in the dynamical structure factor gives rise to a new one- 
dimensional spin liquid state: the so-called non-Haldane spin liquid |(|. 

Clearly, the general classification of all possible one-dimensional spin gapped phases is still 
lacking. However, the integrability of some special field theoretical models has given a lot of 
insight and has been employed to extract the degeneracies, quantum numbers of low-energy 
excited states of spin liquid phases and to provide sometimes the computation of the leading 
asymptotics of correlation functions. This program has been realised for the sine-Gordon model 
associated with the low-energy description of the S=l/2 Heisenberg chain with alternating 
exchange || and the Cu benzoate fl9|,[lO| and also for the O(N) (N=4,6,8) Gross-Neveu model 
for the description of a special frustrated two-leg ladder pT| , weakly two-leg Hubbard ladder 
]f~2| |f~4j , (N,N) armchair Carbon nanotubes and spin-orbital model 15 1. 

The situation is totally different for critical one-dimensional quantum spin systems like the 
spin-1/2 Heisenberg chain where Conformal Field Theory (CFT) allows a complete classification 
of one-dimensional quantum systems or finite temperature two-dimensional classical systems 
jl6|,|l7j . The low-energy spectrum of the lattice model with a continuous symmetry is described 
in terms of representations of a certain current algebra [fl8[ . This affine symmetry determines the 
operator content of the theory and all possible scaling dimensions of the operators are in turn 
fixed by the conformal invariance of the underlying Wess-Zumino-Novikov-Witten (WZNW) 
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model built from the currents of the affine symmetry with the Sugawara construction |lq| . The 
knowledge of the scaling dimensions allows the computation of correlation functions and the 
scaling behaviour of the energies of the ground-state and low- lying states as a function of the size 
of the system pj[ |. The different critical phases are labelled, in particular, by the central charge 
(c) of the underlying Virasoro algebra of the WZNW CFT which fixes the low-temperature 
behaviour of the specific heat of the model . 

A simple example of this description is the Luttinger liquids, a general class of one-dimensional 
interacting electron models whose infrared (IR) behaviour is governed by the Luttinger model 
PH , which have an U(l) affine Kac Moody (KM) symmetry and a central charge c=l. The 
Luttinger liquids are usually described in a bosonized form with the introduction of one gapless 
bosonic field (c=l CFT) and are characterized by single-particle correlation functions having 
branch cuts with non- universal exponents p3]. Such state provides an example of non- Fermi 
liquid behaviour p2[ . The critical theory of quantum spin systems with non-Abelian continuous 
symmetries can also be analysed on the basis of their SU(2)fc affine KM algebra p3| , p4|| . Exam- 
ples of such representation are the critical theory of half-integer spin antifcrromagnets described 
by the SU(2)i WZNW model whereas the family of integrable Hamiltonian of arbitrary spin S 
|p5| belongs to the SU(2)2S universality class This approach has also been employed to 

describe the low-energy physics of Heisenberg antifcrromagnets with a larger symmetry group 
like a SU(N) symmetry |§ on the basis of the SU(N)i WZNW model ||J2^|2§. 

In a recent paper, Andrei, Douglas, and Jerez J30| identified a new non-Fermi-liquid class 
of fixed points, the so-called chiral spin fluids, describing the IR behaviour of one-dimensional 
interacting chiral fermions. The chiral spin liquid (CSL) behaviour appears as a stable critical 
state of one-dimensional interacting fermions with unequal numbers of right-moving and left- 
moving particles. The phenomenological Hamiltonian, in the continuum limit, corresponding 
to this state reads: 




'Hcsl = -iv F > W ar d x R ar - > L J al d x L a i + gJ R ■ J L (1) 



where R) ar (respectively L L) is the creation operator of a right-moving (respectively left- 
moving) fermion with spin index a =t;J. and flavor index r (respectively I): r = 1, 
(respectively I = 1, .., /i). The interaction affects only the spin degrees of freedom and consists 
of a marginal relevant (g > 0) current-current term where the spin currents write in terms of 
the fermions as follows: 

^ }r 

Jr. = 2 Rar a %pRpr 
r=l 

1 H 

ji = - Y. l Ii<^ ( 2 ) 

<r a (a=x,y,z) being the Pauli matrices and (respectively J£) belongs to the SU(2)/ R (re- 
spectively S\J(2)f L ) KM algebra. The Hamiltonian (|l|), in the spin sector, describes thus two 
SU(2) f R and SU(2) f L WZNW models marginally coupled by a current-current interaction. For 
Ir 7^ Jl, i.e. breakdown of the R — > L symmetry implying breakdown of the time-reversal 
symmetry, the Hamiltonian (Q) belongs to a new class of problems not encountered in the 
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above mentionned problems of one-dimensional physics. By a combination of CFT argument 
and Bethe Ansatz techniques, the authors of Ref. |3(]] showed that the Hamiltonian (|l|) for 
fa 7^ II flows in the IR limit to an intermediate fixed point (chiral fixed point) characterizing 
by universal exponents in the leading asymptotic of the electronic Green's functions. From 
the electronic point of view, such a state provides a non-trivial example of a non-Fermi-liquid 
state which differs from the Luttinger liquids: the CSL state. A typical example of possible 
realizations of CSL, given in Ref. |3(|, is a contact between the edges of two integer Quantum 
Hall Effect systems with different filing factors: vl = II vr = /r. 

In this paper, we shall consider another route where the total Hamiltonian (Ti) of a possible 
realization of CSL is time-reversal invariant. Indeed, in the continuum limit, the leading part of 
the Hamiltonian that imposes the strong coupling behaviour decomposes into two commuting 
marginal chirally asymmetric parts H.\ and Ti.2'. 

H = H!+H 2} [Hi,H 2 } = 0, (3) 

where Hi and Ti.2 are not separately time-reversal invariant but they transform into each other 
under the symmetry t — > —t. The model (||) as a whole will display a CSL behaviour in the 
IR limit. Such scenario suggests that other realizations of this state can be found in ladder 
problems. However, in the continuum limit of spin ladders, one has, apart from marginal 
current-current interactions, backscattering terms n a n b (n a being the staggered magnetizations 
of a chain of index a) which have scaling dimension I at the ultraviolet (UV) fixed point and thus 
represent strongly relevant perturbations. Therefore, one has clearly to kill this contribution 
to have a marginal perturbation and thus to find some specific realizations of the CSL state. 
We shall follow, in this paper, two different routes to find some CSL behaviour in the IR 
limit of ladder problems. On the one hand, frustration may play its role either by suppressing 
geometrically the backscattering contribution as in zigzag ladders or by allowing several 

independent coupling constants and the possibility to eliminate the unwanted term by a fine 
tuning mechanism |3^|nj]. In that case, as we shall see later, the CSL found at this special 
point will be destabilized upon switching on weak backscattering contributions but there will 
exist an intermediate low-energy region governed by the CSL behaviour. On the other hand, 
doping effects can also be useful to suppress the backscattering term by some incommensurate 
effects: In the continuum limit, the unwanted contribution will become a strongly oscillating 
perturbation that can be ignored in the long distance limit. 

The remainder of the paper is organized as follows: Different possible realizations of the 
CSL state are introduced in Section II: the (N+l)-chain cylinder model, its asymmetric doped 
version, and the N-channel Kondo-Heisenberg model away from half-filling. The nature of the 
IR fixed point (chiral fixed point) for these models is determined in Section III. In Section IV, 
we present a Toulouse point solution in the N = 2 case that captures all universal properties 
of the chiral fixed point. In that case, the massless degrees of freedom are described in terms 
of an effective S=I/2 Heisenberg spin chain and two critical Ising models. A brief summary 
of this approach has already been published We shall provide here and in Section V 

the technical details of the Toulouse point solution and in particular the determination of the 
physical properties of the 3-chain cylinder model. The nature of the phase diagram of the doped 
models in the N=2 case will be addressed in a forthcoming publication {jfjj. In section VI, the 
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stability of the chiral fixed point of the 3-chain cylinder model is analysed under the presence 
of weak backscattering perturbations. Finally, our concluding remarks are presented in Section 
VII. Three Appendixes (A,B,C) give further technical details necessary for the computation of 
the leading asymptotics of the correlation functions occuring in Section V whereas in the last 
one the stability of the Toulouse limit solution is studied. 



II. POSSIBLE REALIZATIONS OF CHIRAL SPIN LIQUID BEHAVIOUR 



In this section, we shall present some lattice models that are described in the continuum 
limit by an Hamiltonian with only marginal current-current interactions. In the next section, 
we shall show that the low-energy physics of this Hamiltonian is governed by the chiral fixed 
point and thus represents some specific realizations of a CSL state. 



A. The (N+l)-chain cylinder model 



The first model that we consider below can be visualized as follows. It consists of N identical 
S=l/2 antiferromagnetic Heisenberg chains on the surface of a cylinder parallel to each other 
and to the axis of the cylinder. An additional S=l/2 Heisenberg chain, that will be called 
the zeroth chain or central chain in the following, resides inside the cylinder along its axis. 
The surface chains do not interact with each other but all of them are coupled to the central 
one by on- rung interchain interaction (Jj_) and an interchain interaction along diagonals of the 
elementary plaquettes (J x ). The Hamiltonian corresponding to this model is given by: 

N 

H = J2K + H int , (4) 

a=0 

where TiP a is the Hamiltonian of a S=l/2 antiferromagnetic Heisenberg spin chain of index a: 

Hl = Jh^Soj-Soj+i, J,, >0 (5) 

j 

S tt j being S=l/2 spin operators on each site j of the chain of index a = 0,1,.., A. These 
Heisenberg chains are coupled with the interaction: 

N N 

Hint = Jx ^2 ^2 ' + ^ x i^a.j ■ So,j + l + S Q j + i • So,j] ■ (6) 

j a=l j a=l 

In the following, it will be assumed that all interchain exchange interactions are positive and 
weak: < J±, J x <C J\\ to describe the model by a continuum limit approach. Due to the 
presence of the plaquette interaction J x , the cylinder model belongs to the class of frustrated 
spin ladders. For N = 1, this model has been investigated by several groups |37],[ll| and reveals 
no CSL state. As it will be discussed in Section III, the (N+l)-chain cylinder model with N > 1 
is a better candidate for displaying non-trivial CSL physics. 

Let us first consider the continuum limit of the (N+l)-chain cylinder model. The continuum 
description of the S=l/2 antiferromagnetic Heisenberg spin chain is based on the Sugawara 
representation of the SU(2)i WZNW model with a marginal irrelevant perturbation ]23|j24j] : 
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^a = ^ (Jafl ' JaR. + JaL ■ J«l) + l^aR • JaL, a = 0, 1, .., N (7) 

with 7 < and 3 a R,L are the right and left spin currents belonging to the SU(2)i KM algebra 



|l8| , |l7| , p8 39 1 and verify the following operator product expansion (OPE): 

JaL («) JbL H ~ TT-T —.JaL H + — rj 

27T (Z - W) 87T 2 (z - W) 

J afl (^) 4fl M ~ -A H + 



2tt(2-u;) aKV 8^ 2 (*-«;) 



2 



(*) M ~ 0, a,p = x, y, z, a,b = 0, 1, ..,N (8) 



where z = vt + ix, r being the imaginary time, and v ~ Jyao denotes the spin velocity on each 
chain (ao being the lattice spacing). 

The next step to achieve the continuum limit of the Hamiltonian (^J) is to use the continuum 
representation of the spin densities of each chain gj] : 

^£ - S a (z) = J a (x) + (-l) x/ao n a (*) , a = 0, 1, .., N (9) 

with a; = j'ciq and J a , n a denote respectively the uniform and staggered part of the spin density. 
The smooth part J a can be expressed as: 



>aR + JaL, 



(10) 



whereas the n field is the vector part of the primary field of the SU(2)i WZNW model trans- 
forming according to the fundamental representation of SU(2). Assuming weak interchain 
couplings Jj_, J x -C J\\ , one can perform the continuum limit of the cylinder model using Eq. 
(0) and the continuum description @ to obtain: 

H = ^ ( J afi ' Ja/? + JaL ' JaL) + 7 ^ J aR ' JaL + ^ n • n + + flrJ • I, (11) 

a=0 a=0 

where we have discarded all irrelevant operators and oscillatory terms. In Eq. ([ll]), I (re- 
spectively n_|_) is the sum of all uniform (respectively staggered) magnetizations of the surface 
chains: 

N 



a=l 
N 

n+ = ^n a . (12) 



a=l 



Moreover, the bare coupling constants in (11) are given by: 

g = a Q (Ji + 2J X ) 

g = a {Jx-2J x ). (13) 

The interacting terms in ( |ll| ) are of different nature. A first one with the coupling constant 
g is a relevant perturbation of scaling dimension d = 1 and corresponds to the usual interchain 
interaction of non-frustrated spin ladders. A second one with the coupling g describes an 
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interaction between the total spin current of the surface chains and that of the central chain. 
This interaction is marginal relevant and, as long as g is not too small, can be discarded. As 
a result, for generic values of g and g, the low-energy physics of the model will be essentially 
governed by the backscattering contribution. In that case, frustration will play no role except 
for renormalization of spin velocities and mass gaps of the modes that will become massive 
in the IR limit. However, it is important to stress that, in constrast with non-frustrated spin 
ladders, the two coupling constants g,g can vary independently in the model, and there exists 
a vicinity of the line J± — 2 J x (g = 0) where frustration shows up in a nontrivial way. Along 
this line, the low-energy properties of the model are mainly determined by current-current 
interactions: 

a=0 o=0 

which displays in the IR limit CSL physics as it will be shown in Section III. 

B. The (N+l)-chain cylinder model with doped central chain 

As already stated in the Introduction, the effect of doping might be useful to suppress the 
backscattering perturbation in the long distance limit. In particular, one can, for example, 
consider the previous (N+l)-chain cylinder model where some concentration of holes are in- 
troduced in the central chain. Clearly, the very assumption that there can exist such model 
with only one chain doped needs justification. Probably, the simplest one would be to assume 
that the intrasite single-electron energies on the central (eo) and surface (eg) chains are much 
larger than the interchain hopping amplitude: ij_ <C |eo — £s\ so that different chains are out 
of resonance. Anyway, we shall assume in the following that one can dope the central chain 
only. In that case, this special model belongs to the class of one-dimensional electron liquid in 
an antiferromagnetic environment [p0|— p3|] made by the surface chains. 

Due to the presence of doping, the continuum description of the spin density of the central 
is no longer given by Eq. (|^) but takes now the following form (see for instance chapter 16 of 
Ref. §§): 

S (x) = J (x) + cos (2k Ffi x + V27r$ c (xfj n (x) , (15) 

where $ c is a bosonic field associated with the charge degrees of freedom, kp^ being the 
Fermi momentum of the central chain. In the continuum limit, the backscattering contribution 
(no • n + ) acquires now an oscillating factor ~ exp(i(2/cp.o — ir)x) and thus becomes suppressed 
in the long distance limit if the concentration of holes is sufficiently high. As a result, the 
Hamiltonian of this asymmetric doped model in the spin sector is given by Eq. ( |l4"| ) with only 
current-current interactions. 

C. The multichannel Kondo-Heisenberg model away from half-filling 

The last example that we shall give in this paper is the so-called Kondo-Heisenberg model 
with N-channel. This model consists of N identical Hubbard chains (Ti-u) interacting via a 



G 



Kondo coupling with a periodic array of localized spins described by the spin-1/2 operator So,i 
|}l],fl4| [h|. The Hamiltonian of this system reads as follows: 

H = Hu + Jk ^ So,i • S c ,i + Jh ^ S ,i • S ,,+i (16) 



with: 



n 



U 

where c na _i denotes the electron annihilation operator on the site i, and a =t,|, n = 1,..,N 
are respectively the electron spin and channel indices. The conduction electron spin operator 
S C) j at the ith site is defined by: 

S h = 2 Y clia^apCnp^, « = X, J/, Z. (18) 

In Eq. ([l6|), the interaction between the electronic degrees of freedom and the localized spin 
is made by a Kondo coupling with the constant Jk > 0. The Hamiltonian (jl^) can also 
be viewed as a generalization of the multichannel Kondo lattice with an exchange interaction 
Jh > between the impurities spins (simple RKKY interaction). 

In the weak coupling limit when Jk,U <C t,Jn, this Hamiltonian can be represented in 
a bosonized form following the same route as in Refs. (3^^^] for N = 1. The continuum 
description of the localized spin operator (So,j) is still given by Eq. (J9J) whereas that of the 
conduction spin density operator can be expressed in terms of N bosonic fields <I> ac and N 
SU(2)i spin currents 3 a R,L ■ 

N 

S c (X) = (Jafl (X) + J aL (x) + COS (2k F X + V^^ac (a)) n Q (x)) . (19) 

a=l 

With these relations at hand, one can derive the continuum limit of the Hamiltonian jl^ ) 
away from half-filling (incommensurate filling). In that case, the 2kp oscillation of the spin 
conduction operator becomes incommensurate with the alternating localized spin operator so 
that far away from half-filling, only marginal current-current interaction survives in the Kondo 
interaction in the long distance limit: 

2nv N 2irv 

T~L = g " ^ (Jail ' Jail + JaL • Jai) H ^~ (Jo_R ' Jo_R + JoL ' J()l) 

a— 1 

JV 



Jafl ' JaL + 7 Joii • Joi (20) 



a=l 



where I = J2a=i v s ~ 2ta - Ua /2ir, v ~ Jfja , g ~ Jxa , and 7,7' < {U,J H > 
0). Therefore, the multichannel version of the Kondo- Heisenberg model away from half-filling 
belongs to the class of models with marginally coupled current-current interactions. One should 
note that it is more transparent to use a different basis, a spin-charge-flavor decomposition 
^H,^4|, to single out the charge degrees of freedom. In that case, the continuum description 
of the Hamiltonian of the N-channel Kondo-Heisenberg model away from half filling reads as 
follows in the spin sector: 
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T~L S — jy ^ Q R ' ^ R + ^ L ' ^ 3^ ^ 0R ' ^ 0R + ^ 0L ' "J° L ) 



-gj ■ I + 7l fl • II + 7 Jofl • Jol (21) 



III. IDENTIFICATION OF THE INFRARED FIXED POINT 

We investigate now the IR properties of the Hamiltonian ([l4]) with only current-current in- 
teractions which represents the continuum limit of several lattice models as seen in the previous 
section. In particular, the nature of the IR fixed point that governed the low-energy physics 
of these models for N > 1 will be identified and belongs to the universality class of chirally 
stabilized fluids introduced in Ref. |30|| . 

Let us first look at the problem from a perturbative point of view near the UV fixed point 
where all chains are decoupled so that the model has the [SU(2)i]^ +1 ® [SU(2)i]^ +1 symmetry. 
It is conformally invariant with a total central charge: cjjv = N + 1 (N + 1 gapless bosonic 
modes). The one- loop Renormalization Group (RG) equations for the coupling constants g and 
7 are given by: 

g?7 7 2 

v 

(22) 



d In L 2ttv 

dg g 2 



d In L 2ttv 

The bare coupling constant 7 is negative in all the models presented in Section II, so that the 
perturbation 7 ^ a 3 a ^ ■ J a £ is marginally irrelevant and can therefore be neglected: it leads to 
logarithmic corrections. On the other hand, the effective interchain coupling g increases upon 
renormalization {g > 0). Usually, the development of a strong coupling regime (g(L) — > +00) 
is accompanied by a dynamical mass generation and the loss of conformal invariance as for 
the non-linear sigma model or the Gross-Neveu model in 1+1 dimensions. This is indeed the 
case for N — 1 which represents a continuum version of the two-leg zigzag ladder |3l],|3^] in the 
absence of the so-called twist perturbation [j33f . However, for N > 2, the interaction flows to an 
intermediate fixed point (<?* < 00) where the system displays critical properties with a smaller 
(and generally non-integer) central charge than cuv according to the Zamolodchikov c-theorem 
& 

To this end, let us rewrite the Hamiltonian (|l4|) in terms of the leading part which governs 
the strong coupling behaviour: 

27rwi N 2ttv 

7~L = — ^ (Jafl • Jafl + JaL ' JaL) H jj — (Jo_R ' Joii + Joi 1 Jol) + .9 (^OB. ' li + JoL 1 Ir) 

a— 1 

(23) 

where we have omitted the marginal irrelevant contribution and the interaction between the 
currents of the same chirality whose effect will be effectively taken into account by allowing 
the surface chains velocity (v\) to be different from that of the central chain (vo). We shall 
cxplicitely check for N — 2, in Appendix D, by our Toulouse point approach that all these 
discarded interactions will not destabilize the intermediate fixed point. 
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The structure (|23| ) for N > 1 suggests that some degrees of freedom do not participate 
to the interaction and thus will decouple and remain massless. Indeed, the free part of the 
Hamiltonian ( p3| ) is a sum of (N + 1) critical SU(2)i WZNW models whereas the interaction 
describes a coupling between the SU(2)i spin current of the central chain (Jqr.l) and the total 
surface current (J.r,l) which is a SU(2)jv current being the sum of N SU(2)i currents (see Eq. 
(|l2])). To identify the nature of the non-interacting degrees of freedom, we use the following 
decomposition: 

N 

Y[SU(2)i\i = SU(2) N ®g N (24) 
i=l 

where Gn is some piece. The central charge of the critical model with symmetry Gn is: 

c GN = N - c su(2)N =N- — = 4^ (25) 

where we have used that the central charge of the SV(2)k WZNW model is: c = 3k/(k + 2) 
Jl8[ |. Though it requires a proof, for N — 2m, Gn can be viewed as a direct product of m Z 2m - 
symmetric critical models E§] whereas for N = 2m + 1 , it is the product of m Z2,„-symmetric 
critical models and the minimal model A4n+i- We recall here that the central charge of the 
minimal model series M. p with p > 2 is: c p = 1 — 6/p(p + 1) [^6|Jl7| . The previous identification 
reproduces the value of the central charge and also all primary fields have the correct scaling 
dimensions but it is not sufficient to prove the previous CFT embedding. Anyway, in this work, 
we only need the actual value of the central charge (|25|) of the critical model with symmetry 
Gn which one can note that it coincides with the sum of the central charge of the N — 1 first 
minimal models: 



N (N - 1) v^ 1 / 6 



eg* 



y i - —7^ — \ • (26) 

iV + 2 ^ V m(m + l)J v ' 

rn=2 v \ ' 1/ 



The actual model under consideration simplifies thus in the following way: 

H = Hg N +H (27) 

where Tig N is the Hamiltonian corresponding to the critical degrees of freedom with symmetry 
Gn that decouple from the interaction. All non-trivial physics is incorporated in the current 
dependent part of the Hamiltonian (7Y) which can be decomposed into a sum of two commuting 
and chirally asymmetric parts: 

H = H 1 +H 2 , ([Wi,W a ] = 0) (28) 

where 

_ 27TO1 27TOQ 

~~ N + 2 3~ ' ,<? ' 

H 2 = i L . I L + ^ J QR . J 0R + g l L . J 0R . (29) 

The Hamiltonian Hi thus describes two marginally coupled [SU(2)jv]^ and [SU(2)i]l WZNW 
models and TL2 is obtained from Ti\ by interchanging the indices R and L. 
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The properties of chirally asymmetric models with the structure of 7i\ in Eq. (|29|) but with a 
more general symmetry group [SU(2)/ r ].r ® [SU(2)/ L ]i (Jr > Jl) have been studied by Andrei 
et al. |3(J (see Eq. (jl])). These authors argued that such models flow in the IR limit to an 
intermediate fixed point g* (chiral fixed point) whose symmetry is determined by a WZNW 
coset II: 



[SU(2) fR . fL } 1 



SU(2) fL x SU(2) fH - fL 
SU(2) fR 



(30) 



This fixed point has been found by the authors of Ref. |3Cj using the fact that the differences 
between the chiral components of the central charge (cr — c£) and the KM levels (/r — fa) are 
preserved under the RG flow. Morevover, since the total central charge (cr + cl) always de- 
creases upon renormalization according to the Zamolodchikov c-theorem ^7|] , the identification 
of the IR fixed point is equivalent to find a conformally invariant model with the lowest total 
central charge consistent with the fixed values of cr — Cl and Jr — fa. It turns out that the 
solution for this problem is given by the coset model (|30|). The structure of this fixed point has 
also been checked by the computation of the low temperature behaviour of the specific heat 
from the Thermodynamical Bethe ansatz approach (3^] since this physical quantity is a direct 
probe of the value of the central charge of the IR fixed point |^0| . 

Using the result ( |30| ) , one can immediately deduce the structure of the IR fixed point of the 
Hamiltonian (128): 



[SU(2) 



N-1\r 



SU(2) 1 x SU(2) N . 



SU(2) 



N 



[SU(2) 



N-1\ L 



SU(2) 1 x SU(2) N . 



SU{2) 



N 



(31) 



where the first (respectively second) part is the symmetry of the IR fixed point associated with 
Til (respectively H.2)- The central charge in the IR limit (cir) corresponding to the model 
described by the Hamiltonian (fl) of Eq. (E8ft is thus given by: 



cm (N) 



2 (N - 1) (27V + 5) 



(32) 



(JV + l)(JV + 2) 

From this result, we deduce the value of the total central at the IR fixed point of the original 
Hamiltonian (|2^) taking into account the contribution of the massless degrees of freedom (EH) 
that do not participate to the interaction: 

(N — 1) (TV 2 + 5iV + 10) 



CIR (N) = C IR (N) + Cg, 



(33) 



(JV + l)(JV + 2) 

In summary, the low-energy physics of the (N+l)-chain cylinder model, its asymmetric doped 
version, and the multichannel Kondo-Heisenberg model away from half-filling, are governed by 
an intermediate IR fixed point displaying CSL behaviour with a central charge which is not 
integer but rational for a generic value of N. The first two models, in the spin sector, have a 
central charge given by Eq. (^3|) whereas the gapless modes corresponding to the spin sector 
of the N-channel Kondo-Heisenberg model are fixed by Eq. (|3^). Moreover, from Eqs. ( |32| , j33"l ) , 
one can notice that the central charges vanishe for N = 1 as it should be (spin gap phase) 
whereas for N = 2 one has the simple values: cjr(2) = 3/2,0/^(2) = 2 which means that a 
simple and independent approach should be possible in that case. 
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In the following, we shall focus on this N = 2 special case by presenting an exact solution 
of the U(l) version of the N — 2 problem using a Toulouse point solution. This approach has 
been extremely fruitful especially in Quantum Impurity problems and gives a simple description 
of the two-channel Kondo model and also of the Kondo lattice [50|-|5^]. In particular, for the 
two-channel Kondo problem, Emery and Kivelson |50| identified the residual zero point entropy 
stemming from the decoupling of a Majorana fermion degrees of freedom using a mapping 
of the model onto free fermion theory for a very special value of the interaction (Toulouse 
point). Although the position of the solvable point is non- universal, the Toulouse limit solution 
captures the physical and universal properties of the low-temperature behaviour of the model. 
In our problem, the Toulouse point approach provides a direct and transparent reading of the 
spectrum of the model for N = 2. It also gives a non-perturbative basis from which all leading 
asymptotics of spin-spin correlations for the 3-chain cylinder model can be determined. 



IV. THE TOULOUSE POINT SOLUTION IN THE N=2 CASE 

In this section, we present the solution, by a Toulouse point approach, of the model (^) with 
N = 2: 

7~C = g 1 (Jlfi-Jlfi + JlL-JlL + ^2R-^2R + ^2L-^2l) H ^— {3qr.3qr + JoL-Joi) 

+ g (I^.JoL+lL.Jcm). (34) 



A. Decoupling of a Z2 non-magnetic excitation 

As discussed in section III A, the structure of the interaction of Eq. ([m|) suggests an alterna- 
tive representation of the two SU(2)i WZNW models of the surface chains in a SU(2)(g> Z 2 way 
reflecting the global spin rotational symmetry as well as the discrete Z2 related to the inter- 
change symmetry (1 — ► 2). The critical model with the symmetry group Q2 in the description 
( |24| ) for N — 2 identifies with a single critical Ising model (or a free Majorana fermion in the 
continuum limit) since on has from the Goddard-Kent-Olive construction j4Sj ] : 

SU{2) 1 ®SU{2)i ~SU(2) 2 ®Z 2 . (35) 

A simple way to take into account this non-magnetic Z2 excitation (Ising degrees of freedom) 
is to use the representation of two SU(2)i spin currents in terms of four Majorana fermions £° 
and f J53|3|g: 



Iq — Jia + J2Q — — — A |o, a — R,L, 

3 la - ^2a = i faCai " = L (36) 

where the four right and left moving Majorana fermions have the following anticommutation 
relations (a, b = 0, 1, 2, 3): 

{&(*), &(V)}=0 (37) 
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and are normalized according to: 



27T (z — W) 



27T (z — U>) 

S(^)dH~0. (38) 

Using the correspondence ([36|), the Hamiltonian ( |34] ) simplifies as 

H = -i| ^ + Jo] (39) 

where the first part is the free Hamiltonian of the Majorana fermion £° with central charge 
c = 1/2 and 

n = ^Y (Ifl-Ifl + Il-Il) + (Jofl-Jofl + Jol-Jol) + .9 (I«.Jo£ + Ii-Jofl) ■ (40) 

All non-trivial physics is incorporated in the current dependent part of the Hamiltonian (7i) 
which decomposes into a sum of two commuting and chirally asymmetric parts as in Eqs. 
@) for N = 2: 

n^n 1 + n 2 , ([h 1 ,h 2 ] = o) (4i) 

-U - T T I 2nV V T T . T T 

111 — *-R ■ lR H g — JOL ' JOL + g lR ■ JOL 

W 2 = ^j- I L ■ I L + J fl • Joii + 5 II • Jofl- (42) 



B. Identification of the Toulouse point 

We consider now the SU(2) broken version of model ( |4l| ) (g — ► <7|i,<j , .j_) where we allow gu to 
be different from gj_: 

Hi = ^ Ir Ir + ^ Jol-Jol + 3|| I'rJSl + Y {iPol + H.c) (43) 

and 7^2 is obtained from TLi by inverting chiralities of all the spin currents. 

In this section, we present an exact solution of the Hamiltonian (^) using Abelian bosoniza- 
tion. The solution is based on a mapping onto Majorana fermions and exploits the existence 
of a Toulouse-like, exactly solvable point, at a special value (though non-universal) of g\\ where 
the fermions are free. 

The scaling equations for the model (^3|) are of the Kosterlitz-Thouless form: 

d 9\\ = g\ dg ± = g n g± 

d\nL 7r (v\ + vo) ' d\nL ir (v\ + vq) ' 

indicating the increase of the coupling constants upon renormalization for gii > — | • Since 
both the exactly solvable point and the isotropic strong-coupling separatrix gu — \g±_\ occur 
within this range, the exact solution of the anisotropic model (f43| ) is expected to exhibit generic 
properties of the original, SU(2) symmetric, model (^). 



12 



The starting point of our approach is the Abelian bosonization of the component of the 
SU(2)i current Jo associated with the spin excitation of the central chain. Wc introduce a 
massless bosonic field if and write (see for instance Appendix A of Ref. J54j or chapter 13 of 
Ref. f§): 

J QR = ~Jk= d x<PR> J 0L = ~7== d xVL 



It is a simple matter to show that the representation ( ^B|) indeed reproduces the OPEs (||) of 
the SU(2)i WZNW model 

On the other hand, the SU(2) 2 current I can be expressed in terms of a triplet of massless 
Majorana fermions £ a , a=l,2,3 (see Eq. (|3^)). Since we consider the U(l)®i?2-symmetric 
version of the model, it is natural to discriminate (C 1 ,^ 2 ) from the £ 3 field. The two Majorana 
fields (C 1 ,? 2 ) can be combined to form a single Dirac field x which can in turn be bosonized 
with the introduction of a massless bosonic field <&: 

= £g + _ K JV^$r 
XR V2 V^o" 

where the anticommutation relation between \R and XL is insured by the choice [<&r, $l] = i/4. 
On the other hand, to take into account the correct anticommutation with the third Majorana 
fermion £ 3 , one needs to introduce an additional real fermionic degree of freedom k (k 2 = 1). 
Using Eqs. (|3(| Hq), the components of the SU(2) 2 spin current I transform to: 

VTT 

Ir = — — Cfl«e Ti ^* R , if = — = ei«e ±iv ^ $i . (47) 

The chiral aymmetric parts of the total Hamiltonian TL ([i"l| ) can then be written in the following 
bosonized form: 

Hx = v {d xVL f + v x (d x $ R ) 2 - iy £ 3 R d x & 

+-?Ld x ip L d x $ R + l l±—£_ RKC os(^<$> R + V^ip L ) (48) 

V2tt 2(7ra ) 1 V 7 



and 



' d x <p R d x $ L + '-i^-^tlKcosfV^QL + V^m)- (49) 



V2tt 2 (7ra ) 

One can eliminate the cross terms d x ipd x $> by performing a canonical transformation: 




cha sha 
sha cha 




(50) 
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and 



PR 



cha sha 
sha cha 



$1L 
®2R 



with 



th2a = 



■911 



Under the transformation (|50| ) and (^3l|), the arguments of the cosines in Eqs. (| 

(VZcha + sha^j $ 2 L + (cha + \/2sha) <&ir 
\/2cha + sha ) § 2 r + (cha + \/2sha ) $il 



47r$ L ^ 



877(^9^ 
877^^ - 



4/r 

47T 



We immediately observe that for a special value of a given by 

tha = — , 

V2 



(51) 
(52) 

19]) become: 
(53) 
(54) 

(55) 



we are at the free fermion point where the two cosine terms acquire respectively conformal 
weights (0, 1/2) and (1/2,0) and can be expressed in terms of the left and right components of 
a single Dirac fermion tp. When a satisfies both Eqs. (52, 5^), g\\ takes a special positive value 
gt (the Toulouse point): 



4tt 



9\\ = "3 

and the total Hamiltonian (it) reads as follows 



(56) 



H 



+ 



2 (7ra ) 



3/2 



£|k cos ( 



47T$ 



2L 



.. 2 i) ! '••'•">( 



U2 

2 



(9^2) 2 + (^e 2 ) 2 



4tt$ 



(57) 



2 (7ra ) 

where we have reestablished the complete bosonic fields i> a = $ q l + (a = 1, 2) and their 
dual counterpart Q = <& a £ — <& a fl ; . The associated velocities u\ and u 2 are given at the Toulouse 
point by: 



Ul 



U2 



(58) 



3 ' 3 

The Toulouse point solution is thus stable provided that all the velocities of the modes are 
positive, i.e. 1/2 < vq/vi < 2. Since we are at the free fermion point, we can refermionizc 
the two cosine terms in Eq. (57). Introducing a pair of Majorana fields rj and £ and using the 
correspondence: 

'//.'./ + Kr.l 



s/2 



;i ±l V / 47r$2H,. 



the two cosine terms become simply: 

K COS ( \/47T$2_R,L J = TjR.L. 



(59) 



(60) 
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Therefore, the Hamiltonian (Ti) takes the following final form: 



n = y [(^*i) 2 + (^©i) 2 ] - 1J y iCnd^R - ( L d x a] - y ~ 
1U2 



,, [mdxVR - VlOxTil] + im [£%t}l - m^l] (61) 

with m = gj_/2a a 7r. Taking into account the contribution of the Z 2 Majorana fermion £° which 
has decoupled from the beginning, the original Hamiltonian ( p^ ) reads at the Toulouse point: 

w = Y [(^i) 2 + (d x ei) 2 ] ~ ^ [Ck&Ch - a&a] - y" - £^S] 
- y" - ci^ci] - ^ - + im - • ( 62 ) 

At the Toulouse point, the massless degrees of freedom are thus described in terms of an 
effective S=l/2 Heisenberg spin chain associated with the $1 field and two decoupled critical 
Ising models Q so that the total central charge in the IR limit is: cm = 2 in full agreement 
with Eq. ( |33| ) for N = 2. One should also notice that the nature of the critical fields found in 
the Toulouse limit approach reproduces, for N = 2, the structure of the symmetry group ( |3~l| ) 
of the fixed point corresponding to the Hamiltonian (Tt). Indeed, substituting the value N = 2 
in Eq. (|3l|), we see using Eq. ([55]) that this symmetry group simplifies as: 

[[SU(2)i] R ® [Z 2 ] L ] <g> [[St/(2)!] L ® [Z 2 ] fl ] (63) 

which is in full agreement with the structure of the Hamiltonian ( |6l| ) at the Toulouse point. 

The remaining part of the Hamiltonian (|6^) has a spectral gap and describe the hybridization 
of the Majorana £ 3 and 77 fields with different chiralities. Indeed, it can be written as a sum of 
two commuting part TL + + 7i_ with 



n- = -^e R d x e R + ^v L d xVL + ^ ~ v L e R ) . (64) 

The resulting spectrum is given by 



E+(k) =k Ul ±\lk 2 [ Vl ~t Vo ) +77i 2 (65) 



g_ (fe) = -km ± ] jk 2 ( ^ 1 ^ Va j +m 2 , (66) 

which, apart from the velocity anisotropy, is reminiscent of the structure of excitations found 
in the Majorana approach for the Kondo lattice |5q| . One should stress that each Majorana 
fermion hybridizes with the other with the opposite chirality and that there is no coupling 
between Majorana fermion of the same chirality. This reflects the chiral nature of the fixed 
point. One should finally note that our solution contains the SU(2) point when g± = g» and 
the mass at the Toulouse point reads in that case: m = gjj /27ra = 2(vi + u )/3ao. 

The exact solution based on the existence of a Toulouse point allows us to get a simple 
understanding of the structure of the critical elementary excitations and to estimate the leading 
asymptotics of the correlation functions. In the following, we shall realize this program by 
directly focusing on the physical properties of the 3-chain cylinder model. 
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V. PHYSICAL PROPERTIES OF THE 3-CHAIN CYLINDER MODEL 



In this section, the spectral properties of the 3-chain cylinder model are discussed and we 
compute the leading asymptotics of all spin-spin correlation functions at the chiral fixed point. 
To this end, we first express the smooth and staggered magnetization of each chain in terms of 
the suitable fields describing the long distance physics at the chiral fixed point. 

A. Representation of the spin currents and staggered spin fields at the Toulouse point 

In this subsection, the spin currents (3 a R,L) and staggered spin fields (n a ) of the chain 
a = 0, 1,2 are related to the two critical Majorana fermions (£°,C) an d the massless bosonic 
field ($1) at the Toulouse point. To get a better insight in the structure of excitations, we 
introduce the smooth part {J a ) and staggered part (N a ) corresponding to the effective S=l/2 
Heisenberg spin chain described by the bosonic <E>i field. Using Eq. (fl5"l), one has the following 
correspondence: 



n 2na 2-Kao 
The staggered part of a SU(2) spin density can also be expressed in terms of the bosonic <I>i 
and its dual field 6i by (see Ref. Q or chapter 13 of Ref. |3S|l): 

A 



Af z = sin (V2^$i) 

nan V / 



7ra 

N ± = — e ±iv ^ §1 (68) 

where A is non-universal constant related to the band gap for the charge excitations of the 
underlying Hubbard model. The WZNW primary operator g transforming according to the 
fundamental representation then writes: 

g= JL=(c + fi.g) (69) 

where e = Acos(\/27r$i)/7rao is the dimerization operator. With the bosonic representation 
( |67| , |68| ), one can check that the primary field g ( |69| ) verifies the following OPE: 

Jl (z) g (u, Q) ~ - _ , ; 1 , . g (u, Q) (70) 



2tt (z - uj) V 2 

J% (z) g (u, CD) ~ - 1 g («, u>) [^) (71) 
2ir (z — to) \ 2 J 

which corresponds (up to the normalization factor l/27r) to the definition of a WZNW primary 
operator of Ref. JlTj] (chapter 15). 



1. Representation of the spin currents 

We begin our analysis by considering the SU(2)i spin currents Jo_r.l of the central chain. 
Using the canonical transformation ( |50| , |5l| ) with the special value (^), we obtain: 
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<p L = V2$ 2i - (72) 



from which together with Eq. ([15|) we deduce the transformation of the spin currents of the 
central chain at the Toulouse point: 

J 0R.L = ~J£,R + —7=dx®2R,L 

T± _ rr± =Fiv / 167r4>2R 

J 0R ~ J L e 

j ol = j£f£ iy/m *™' (73) 

where we have used ( |67j ) to express the spin currents of the central chain in terms of the 
"physical" current J a . Using the correspondence (|5|), d x &2R,L can be related to the Majorana 
fermions T)R t j, and (r,l'- 

d x $ 2 R,L = ^VR,Ur,l- (74) 

Consequently, at the Toulouse point, we deduce the expression of the z-part of the spin current 
of the central chain: 

JSr = + ' 1t ]r(r 

JSl = -Jr + ^lCl. (75) 



Since the bilinear fields ^^lCr.l are short-ranged, we can simplify further ( |75| ) and obtain 
the following correspondence valid for the estimation of the leading contribution of correlation 
functions: 

Jm,L ~ -JLr- (76) 

Notice that the chiralities of the central chain and the physical spin current are opposite. We 
shall return to this important point when discussing the nature of elementary excitations of 
the model. For the transverse part, the bosonic operator exp(±iVl67r$2fl,i) has a non-zero 
expectation value, one thus obtains ignoring short-ranged pieces: 

J§R ~ Ir^L 

j or ~ 1r3l 
Jol ~ 7l J£ 

Jol ~ HJr (77) 



where ^r,l are defined as: 



{j V ^*™- L )=lR,L- (78) 



The actual value of Jr,l is computed in Appendix A where one can show that Jr — Jl = lm 
and 7 m is real. One thus concludes that: 

Jqr,l ~ "<™-Jl,r- (79) 
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The representation of the SU(2)i spin currents of the surface chains (J a R,L, a = 1,2) proceeds 
in the same way. To each chain a, we associate a massless bosonic field tp a and its dual fields 
# a to write the currents as in Eq. (p>|): 



T z — 

J aR — 



1 



S2tt 
1 



dxVaR , 



T z — 



1 



J " R ~ 2.V ( „, 



1 =Flv / 8ir^ a jj t± _ 



/2tt 



27ran 



C 



±lV87T¥'a 



(80) 



To relate these fields to the bosonic field <E> of Eq. ((47|), let us introduce the symmetric and 
antisymmetric combinations: 



¥>± = 4= fax ± Va) , *± = 4= 0? i ± 



(81) 



so that 



1 



%R,L 



2^ 
1 



(d x ip +RtL + T a d x ip- Ri i) 



J « R ~ 27T,/., 



J;' 



1 



ai 2vra " 

where ri = 1 and T2 = — 1. Moreover, since we have the relation: 



(82) 



(83) 



and the expressions ( f47| ) for the bosonization of the SU(2)2 current I, we obtain the following 
correspondence: 



¥>+R,L = ®R,L 



(84) 



The difference between the SU(2)i currents of the surface chains can also be expressed in terms 
of the Majorana fermions £ a ,a = 0,1,2,3 (see Eq. (j36|)). Using Eqs. 82), one thus finds 
the identification: 

dxV-a = iV^Ca^a ) a = R, L 

sin (V^Trip-R^J — y/nao i£rK 

sin ^V^r^-L^ = — (85) 



With these relations at hand, we can now relate the spin currents of the surface chains 
a function of the Majorana fermions £°, £ 3 and the bosonic field $: 

j 1r,l = 2^ dx ® R ' L + 2 Ta ^ R ' L ^ R ' L 



as 



J aR — 



1 



J, 



1 



^±lv / 47r* £ 



'i&K ± T a f R K) 

(i^«±r a ^«). 



(86) 
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Under the canonical transformation (|5(J, pl|) , the Majorana fermions £°,£ 3 remain unchanged 
while the chiral bosonic fields transform at the Toulouse point according to: 

$ L = \/2$ 1L - $ 2i? . (87) 

Therefore, using the relation (|74| ) between the Majorana fields (, r/ and the bosonic field $2, we 
obtain for the z-component: 

JaR =3r~\ {VLCL ~ T a &&) 

and ignoring the short-ranged contributions, we get: 

JaR,L ~ 3r,l- 

In the same way, the transverse parts of the spin currents of the surface chains simplify as 
follows: 

JaR = -™ j£ (mf R ± CUli ± T ar] Ul - ITaCUl) 

JaL = -Vaoj£ (imfL ± Crfl ± Tamil ~ ^G^i) • (90) 

Majorana bilinears, in this equation, built from £ R L and ^ R L are short-ranged and from the 
hybridization (p4|), the only non-zero vacuum expectation values are: 

1r,l = mao{t RL r) L!R ). (91) 

The value of this expectation value is estimated in Appendix A and one finds: j R = Jl — 7m 
(7™ being real). With the same accuracy than in Eq. ([79]), the transverse spin currents of the 
surface chains are thus given by: 

J aR ~ Jr (7m + ia T a 7rCi^) 

JaL ~ Jl (7m + ia Q T a 7rCR£l) ■ (92) 



2. Representation of the staggered spin fields 

Let us first analyse the staggered spin field of the central chain no. As in Eq. (|68|), this field 
can be expressed in terms of the bosonic field ip and its dual field 

A 



sin (\/2tt'p) 

nan V / 



nan 

n ± = A e ±i^. (93) 
na Q 

Using the transfomation of the chiral bosonic fields (|f^), one finds at the Toulouse point: 

if = V2i> 2 - $i 

tf=\/2e 2 + 6i. (94) 
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We deduce then the following result: 

71q = — (sin (V47ri> 2 ) cos (\/27r$i^ - sin (\/27ri>i^ cos (%/47ri> 2 ^ 

+ = Aeiiv^e, e± iv^e 2 _ (95) 
ira Q 

One can simplify these expressions using the correspondence ( |59| ) to obtain the relations: 
cos {y^n^-i^j = iTrao fe^L + (r(l) 
sin (v / 47r$2) = iTrao (~VrCl + CrVl) 

from which, we deduce the leading part of the representation of the staggered spin fields no 
associated with the central chain at the Toulouse point: 

Uq ~ -iira ( R ( L Af z 



n 



± ~ iTTOoCflO^ (97) 



■-^-sin (V^r^a) 
nan V / 



where A/" corresponds to the staggered magnetization of the effective S=l/2 Heisenberg spin 
chain (see Eq. (|68|)). 

The same correspondence for the staggered magnetizations (n a ) of the surface chains is more 
difficult to obtain. The reason for this is that it involves nonlocal strings of the Majorana 
fermions £ a , r\ and (. As we shall see, the trick will be to relate all these quantities to the order 
and disorder parameters of the underlying Ising models present in the problem. The Appendix 
B gives a review of the basic definitions for the bosonization of two Ising models and introduces 
also the different fields which will be useful in this subsection. 

The staggered part (n a , a = 1,2) of the spin density of the surface chains can be expressed 
in terms of the bosonic field <p a and its dual d a as in Eq. (|9J 

_A_ 

ira 

nt = —e ±iV ^*°. (98) 
7ra 

Introducing the symmetric and antisymmetric combinations d8l| ) of the bosonic fields, we get: 

n z a = — (sin (V^*) cos {\fnip-) + r a cos (Vtt*) sin (V^-)) 

nt = A e ±iV^e e ±iwW- ; (99) 

where we have used the fact that (f+R,L identifies to the bosonic fields Qr,l- According to the 
canonical transformation ( p7| ) at the Toulouse point, the bosonic fields <&, 9 transform as: 

$ = V2$i - $2 

e = \/2ei + e2 (too) 



whereas the fields , are not affected by this transformation. Using Eqs. (BS, B9) of 
the Appendix B, one can then relate the staggered spin fields n a of the surface chains to the 
difffcrent order-disorder operators of the underlying Ising models of the problem: 
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(sin (y2n^^j (/i 5 M4M3Mo + 1^05040-300) 



K = 

ira 

+ cos {y^ii^^j (-0-5(74^3^0 + Ta^ma^ao) 
7ra 

±1^5040-3^0 ± ir a cr 5 ^4/^30-0) . (101) 

The Ising models labelled "3,5" are non-critical due to the hybridization d&4| ) of the Majorana 
fermions rj, £ 3 at the Toulouse point. As shown in Appendix C, the order-disorder operators 
(0-3,5, ^3,5) of the corresponding Ising models are short-ranged and have the following vaccum 
expectation values: 

(^3^5) = (0305) = [i ^ 

M = (^305) = 0. (102) 



Consequently, ignoring short-ranged contribution in Eq. (101), we finally obtain the correspon- 
dence between the staggered magnetization of the surface chains and the critical fields at the 
chiral fixed point: 

n l ~ (/U4MO + T a <J 4 <J ) 

n± ~ fj,JV ± (/X4M0 - TaO 4 o ) . (103) 

Let us summarize our results obtained in this subsection. Using our Toulouse point approach, 
the spin currents and staggered fields of the three chains have been expressed in terms of the 
critical fields at the Toulouse point. Up to short-ranged pieces, the spin currents are given by: 



J 0R,L 


J ~$L,R 




J 0R,L 






J aR,L 






J aR 


J 3r {lm ' 


h ia T a ir( L £° R 


J± r 


- Jl (%. - 


h ia T a ir( R £l 



(104) 

whereas the staggered magnetizations of the different chains write as follows 

nl ma CR.C,LM z 

tiq ~ maoC.RC.L.N' 

n z a ~ (jj\f z (/^Mo + T a a 4 a ) 

n± ~ /jM ± (/^Mo - t q o 4 o ) . (105) 

These results enable us to discuss the spectral properties of the 3-chain cylinder model as well as 
the computation of the leading asymptotics of the spin-spin correlation functions at the chiral 
fixed point as described in the following section. 

B. Spectral properties 

There are two different kinds of massless elementary excitations at the Toulouse point: mag- 
netic excitation described by the field $1, and non-magnetic, singlet, excitations associated 
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with the two Majorana fermions £° and £. These two groups of excitations are decoupled at 
the chiral fixed point. Let us first discuss the nature of the magnetic excitations. 



1. Magnetic excitations 

The magnetic elementary excitations of the system correspond to the spinons of the effective 
S=l/2 Heisenberg chain described by the bosonic field These spinons appear only in pairs 
in all physical states and carry a spin S=l/2. ft is very important to notice that, due to the 
mixing of different degrees of freedom reflected in the canonical transformation (50 5T\), the 



"physical" spinons, i.e. those defined as \/tt/ 2-kinks of the field $! should not be misleadingly 
identified as the spinons of the central chain. To get a better understanding of the structure 
of the spin excitations at the chiral fixed point, let us recall the expressions ([75], |88|) of the 
currents J-f 2 and J§ in terms of the "physical" current J z — (l/^/2n)d x ^i at the Toulouse 
point obtained in the previous subsection: 

J 1R(L) = 3r[L) ~ 2 \VL(R)Cl(R) - €r(L)€r(L)) ' 
J 2R(L) = ^R(L) - 2 (vL{R)Cl{R) + > 

J 0R(L) = -Jl{R) + ' 17 1R(L)Cr(L)- (106) 

At energies |cj| 3> m, where all Majorana fields can be considered as massless, Eqs. ( |10G| ) 
transform back to the standard definitions of the currents of the three decoupled chains. In this 
UV limit, one has a picture of three groups of independently propagating spinons. However, 



in the IR limit Qu>\ <C m), all Majorana bilinears in Eq. (106) are characterized by short- 
ranged correlations (since the Majorana fermions r\ and £ 3 are massive), implying that strongly 
fluctuating parts of the currents of individual chains are no longer independent. The spinons 
of the individual chains turn out to be strongly correlated in the IR limit to participate at the 
formation of a single, physical, current J z . In fact, the physical spinon represents a chirally 
asymmetric, strongly correlated state of three spinons. Indeed, consider, for instance, a right- 
moving -y/7r/2-kink of the field $i, representing a physical spinon with the spin projection 
S z = 1/2. According to the exact relation 

^fl(L) = JlR(L) + ^2R(L) + JoL(R) (107) 

following from Eq. ( |I06| ), such an excitation is a combination of two right-moving spinons of 
the surface chains, each carrying the spin S z = 1/2, and a left-moving antispinon of the central 
chain, with S z = —1/2. The rigidity of such a state is ensured by a finite mass gap in the (77— £ 3 ) 
sector of the model. Therefore, in physical excitations, two spinons of the surface chains form a 
bound state with an antispinon of the central chain. Otherwisely stated, the polarization of the 
long-wavelength magnetic excitations of the surface chains is always the same but opposite to 
the polarization of the central chain. This peculiar structure of the elementary spin excitations 
at the chiral fixed point is clearly reflected by the expression for the spin velocity U\ of the 
bosonic field $1: u\ = {2v\ — «o)/3. 

These excitations correspond to the magnetic excitations of the model since a uniform mag- 
netic field H along the z direction only couples to the $1 field. Indeed, in the continuum limit, 
this magnetic field couples to the smooth part of the total spin density of the three chains: 
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H 

Uh = — 7= (d x ipi + d x (pi + d x Lpi) , (108) 

v27T 



and performing the Toulouse transformation (194, 100), one obtains: 



H H = 7==0 x $i. (109) 

•V/27T 



The uniform susceptibility x of the model corresponds thus to that of a single Hciscnbcrg S=l/2 
spin chain and it is given by (in units of g^s)'- 

13 . , 

Y= = (110) 

2ttu 1 2vr (2vi - v ) y ' 

Notice that due to the non-trivial structure of the spinon at the Toulouse point, the uniform 
susceptibility of each chain does not add in the case of a truly uniform field. In contrast, the 
susceptibilities would add coherently if one considers a staggered magnetic field (H a ) in the 
transverse direction: H a = H(8 a \ + 5 a 2 — S a o), a being the chain index. In that case, one finds 
the contribution: 

xT = (in) 



2. Pseudo-charge degrees of freedom 

Apart from the non-trivial nature of the spinon, the chiral fixed point manifests itself in the 
existence of massless Z 2 singlet excitations. A first non-magnetic excitation stems from the 
Majorana fermion £° that decouples from the rest of the spectrum. This Majorana fermion 
describes collective excitations of singlet pairs formed on the surface chains. The nature of 
the other massless Majorana fermion £ at the chiral fixed point is less transparent: It is a 
highly nonlocal object when expressed in terms of the original spin operators. A simple way to 
understand the role of these singlet excitations is to combine the Majorana £° and ( fields into 
a single Dirac \l/ c fermion and then bosonize it with the introduction of a bosonic field $ c : 

y cR = ± = 1 p iV4^^r 

^ cL = iLyct = _^ =e ^v^i cL (112) 

The corresponding massless bosonic field $ c ressembles the scalar field describing the charge 
degrees of freedom in the one-dimensional SU(2) Hubbard model away from half-filling in the 



limit U — 00. This analogy can be seen by looking at the representation (105) of the staggered 
fields at the Toulouse point. Using the bosonization result of two Ising models (B6), the total 
staggered magnetization n + of the surface chains reads as follows in terms of the bosonic field 

n + = ni + n 2 — 2/z/K4 / Lt A/ 

- cos (V^c) AT. (113) 
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This bosonized version of n + is reminiscent of the staggered magnetization of the Hubbard 
model away from half- filling at U = oo (K c = 1/2 in that case, see for instance Refs. |38| , [39|j ). 
Apart from a velocity anisotropy (vo 7^ 1*2), the fermionic fields £° and C plays the role of the 
charge degrees of freedom in the corresponding Hubbard model. We shall hence refer to the 
gapless Majorana fermions £° and C as "pseudo-charge" excitations which account for the central 
charge c = 1 at the IR fixed point. The spin and pseudo-charge separation is already manifest 
at the Toulouse point in the Hamiltonian ( |62| ) since these fermions fields are decoupled from the 
magnetic excitations described by the bosonic field $1. This separation of the different modes 
will also be clearly seen in the expression of the leading asymptotics of correlations functions 
that we now estimate. 

C. Spin-spin correlation functions at the chiral fixed point 

We are in position to compute the leading asymptotics of the spin-spin correlation functions 
of the 3-chain cylinder model at the chiral fixed point. Let us begin our analysis by considering 
the spin-spin correlation functions between the central spins. 

1. Correlation functions between the central spins 

Since the SU(2) symmetry is broken to an U(l) symmetry in our Toulouse point approach 
(S|| 7^ 9±)i we shall discriminate between the z and perpendicular components of the correlation 
function: 

(5 2 (x, r) S$ (0, 0)) = (JS R (x, t) r 0R (0, 0)) + (R—> L) 

+ (-1)^(^(^)^(0,0)), (114) 

(S+ (x, t) S - (0, 0)> - ( J + (x, t) J m (0, 0)> + (R^L) 

+ (-lT /a °(n+(x,T)n^(0,Q)). (115) 

Using the representation ( |104| , |105| ) of the uniform and staggered parts of the spin density of 
the central chain at the Toulouse point, we get the following leading behaviours: 

(5 z (x, t) (0, 0)) ~ (J% (x, t) J Z R (0, 0)) + (R -> L) 
+ (-lf /ao ^al{N z (x, t) M* (0, 0)) (d (x, t) Cr (0, 0)) (Cl (x, r) £ L (0, 0)) (116) 

and 

(S+ (x, t) So (0, 0)) ~ jl(J+ (x, t) J r (0, 0)) + (R^L) 
+ (-l) x/a " Tr 2 al(N + (x, t) M- (0, 0)) {Cr (x, t) (r (0, 0)) (Cl (x, t) Cl (0, 0)) . (117) 

Accordingly, the estimation of these correlations simply reduces to the computation of two- 
point functions of a massless bosonic field $1 and a massless Majorana fermion C- The final 
expressions for the leading asymptotics of the correlation functions between spins in the central 
chain at the chiral fixed point are thus 
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(s* (x,t)s z (o,o)) ~--L(- — * + 1 



8tt 2 I (x + iuxr) 2 (x — iuit) 2 

2 r 



+ ( 1)X/a " »j (x* + U 2 T 2 ) {X 2 + u 2 t2) 1/2 ' ( 118 ) 



(S+(x,r)Sa(0,0))~-^ 



4tt 2 I (x + iuit) 2 (x — \u\t) 

( _,yx/ao ^ 2q 1 

' \ > a— 2 r„2 i . 



4tT 2 (x 2 +U 2 2 T 2 ) (a;2 +u 2 T 2)l/2 



(119) 



2. Correlation functions between spins of the surface chains 

Let us now compute the spin-spin correlation functions between the surface chains: 
(5* (x, t) Si (0, 0)) = ( J* aR (x, r) J z bR (0, 0)) + (R^L) + (— l) X// °° („» ( X) T ) n \ (0, 0)), (120) 

(5+ r) S fe " (0, 0)) = (J+ (x, r) J bfl (0, 0)) + (R - L) + (-l)^ (n+ (s, r) n b " (0, 0)). (121) 
where a, 6 = (1, 2) is the chain index. We shall follow the same route as for the central chain by 



using the expression of the uniform and staggered spin fields (104, 105) at the Toulouse point 



(S z a (x, r) S§ (0, 0)) ~ {J Z R (x, r) J% (0, 0)) + (JJ - L) 
+ (-l)' T/ao M 2 (AP(x,r)AP(0,0)) 
((/i 4 (x, r) ^ 4 (0, 0)) (mo (a, r) ^ (0, 0)) + T a T b (<j 4 {x, r) a 4 (0, 0)) (<t (a, r) ct (0, 0))) (122) 

where we recall that n = 1,T2 = —1 and the perpendicular part reads 

(5+ (x, r) S 6 " (0, 0)> ~ (J+ (x, r) ^ (0, 0)) (7* + r a r b a 2 (C L (x, r) ( L (0, 0)) (£° (x, r) ^ (0, 0))) 

+ (J2, L) -> (L, R) + (— l) x / a ° /i 2 (AA+ (x, r) AA" (0, 0)) 
(( M 4 (or, r) M4 (0, 0)) (/x (x, r) /x (0, 0)) + r a ^(a 4 (x, r) a 4 (0, 0)) (<r (x, r) <t (0, 0))) . (123) 

The Ising models labelled by "0" and "4" associated with respectively the Majorana fermion 
£° and £ are decoupled and critical at the chiral fixed point. The correlation functions of the 
order and disorder operators 00,4 , ^o,4 & re known exactly. In the long time and long distance 
limit they are given by (see Appendix B): 



(fi (x, r) M o (0, 0)) = (a Q (x, r) a (0, 0)) 

1/4 



(x 2 +w 2 r 2 )i/8' 



(124) 



(^4 (x, r) ^4 (0, 0)) = (0-4 (a:, t) 0-4 (0, 0)) 

1/4 



(X 2 +U 2 2 T 2 )^ 



(125) 



25 



This enables us to obtain the result for the asymptotics of the correlation functions between 
spins of the surface chains: 

<S2 (x, r) SS (0, 0))~~ ( J , 2 + !— rj) 

oTT \(x + lMlT) (x — lUiTj / 



2\2 



+ (-l) X/O0 ^44 



^ 2 V^o (x 2 + u\t 2 ) 1/2 (x 2 + v 2 t 2 ) 1/S (x 2 + u 2 t 2 ) 1/8 ' 



(126) 



(5+(x,r)5 b -(0,0)) 



' 16tt 4 



74 1 , , 1 

47r 2 I (x + iu\T) (x — \u\tY 

111 



(x - hut) (a: - i«ir) (x + iu 2 r) ( x + i Ul r) + i^ir) (x - iu 2 r) 
2^ 2 A 2 1 1 1 



+ (-V X/a °5ab- 



(x 2 + u 2 t 2 ) 1/2 (x 2 + v 2 t 2 ) 1/S (x 2 + ulr 2 ) 1 



(127) 



Here again, the leading contribution of the uniform part of the correlation functions is that 
of a single S=l/2 Heisenberg spin chain. The other critical modes manifest themselves in 
subleading contributions in the uniform part and directly in the staggered part. In particular, 
the contributions of the Majorana fermions £° and C m the staggered part show up through 
their corresponding order and disorder Ising operators. 

3. Correlation functions between spins of the surface chains and of the central chain 



The last spin-spin correlation functions to estimate at the chiral fixed point consists in cor- 
relation between a spin of the surface chain S a (a = 1,2) and a spin of the central chain Sq. 
Using the representation ( |104| , |105[ ) of the different spin fields at the Toulouse point, we obtain 
following the same route as for the previous correlation functions: 



1 



1 



(x + Iuit) 2 (x — iuir) 2 



(128) 



<S+(x,t)S o -(0,0))~- 



7m Jn 

Air 2 



(x + iwir) 2 (x — iuir) 2 



(129) 



It is interesting to notice that the staggered parts of these correlation functions are short-ranged 
at the chiral fixed point. This result stems from the very special structure of the spectrum of 
the Hamiltonian at the Toulouse point discussed above. 

We end this section by estimating the NMR relaxation rate 1/Tj at low temperature. The 
slowest correlation functions of the model at the chiral fixed point correspond to the staggered 



correlations between the spins of the surface chains (see Eqs. (126, 127))- We emphasize 
that the exponents occuring in these correlation functions are universal and characterize a new 



universality class in spin ladders. We deduce from these correlation function (126, |127| ) the 
low-temperature dependence of the NMR relaxation rate: 

1 



=- ~ VT 



(130) 
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in contrast with the S=l/2 Heisenberg spin chain where 1/Ti ~ const |38|39[]. Finally, we 
remark that the result ( |l30| ) corresponds to the leading behaviour of the NMR relaxation 
rate of the one-dimensional Hubbard model away from half-filling at U = 



38 3£ 



The 



main reason for this stems from the presence of the pseudo-charge degrees of freedom £ , £ in 
the correlation function (126, 127) that enter through their corresponding Ising order-disorder 
operators (00,4, /xo,4) with scaling dimension 1/8. 



VI. STABILITY OF THE CHIRAL FIXED POINT OF THE 3-CHAIN CYLINDER 

MODEL 

In Sections IV and V, we have carefully analysed the Hamiltonian consisting in only current- 
current interaction describing two marginally coupled SU(2)2 and SU(2)i WZNW models. The 
model with only current-current interaction has been solved exactly using a Toulouse point 
approach. At the strong coupling fixed point, this Hamiltonian belongs to the chiral stabilized 
liquids universality class with central charge c = 2. The critical fields at the chiral fixed point 
consist of two decoupled gapless modes: a magnetic sector described by an effective S=l/2 
Heisenberg chain associated with the bosonic field $x an d the pseudo-charge degrees of freedom 
stemming from the two Majorana fermions £° and f. Moreover, the Toulouse point solution 
captures all universal properties of the chiral fixed point as shown in Appendix D. 

The next step of our approach is to analyse the stability of this chiral fixed point. The 
non-perturbative basis provided by the Toulouse point approach enables us to investigate the 
stability of the chiral fixed point under several operators such as backscattering contributions. 
In particular, we shall first investigate the effect of the interchain backscattering operator to 
deduce the nature of the low energy physics of the 3-chain cylinder model in the vinicity of the 
line g — 0. 



A. Effect of the interchain backscattering perturbation at the chiral fixed point 

Let us begin by studying the effect of the interchain coupling in the 3-chain cylinder model 
in the vinicity of the chiral fixed point: 

H b =gn -n + . (131) 

From Eq. (|f), we see that the coupling constant g is independent from the coupling g of the 
current-current interaction. One can thus always suppose that g is small to investigate the 



backscattering term (131) as a weak perturbation at the chiral fixed point: \g\ -C 1. Using the 



bosonized form of the staggered spin fields given by Eqs. ( |93| , Pq ), we have: 
H b = ( cos (Vnd-) cos (v 7 ^ - a/tt©) + cos (y/n<p-) sin (%/2^) sin (v^)) (132) 



so that in terms of the different fields at the Toulouse point, the backscattering term (131) 
reads as follows: 
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Kb = -^f~2 { cos (-v/th?-) cos (^©2) - \ cos (y/nip_) COS (^$2) 



+ icos( % /¥ ¥ ?_)cos(3V^>2- %/8^$i)^ ■ (133) 

Using the Ising dictionnary (see Appendix B), it can be shown that the contribution in the 
magnetic sector of Eq. (133) is irrelevant. The backscattering term affects mostly the pseudo- 
charge sector. Indeed, using the results ( |C9| ) from the structure of the massive modes at the 
Toulouse point, one finds the following estimate at the chiral fixed point: 

n b ~ ^/f (134) 

We thus conclude that the backscattering operator opens a gap, A c , in the pseudo-charge sector 
but has no effect on the magnetic (spinons) excitations. Standard scaling arguments give an 
estimate of the mass gap since the operator ( |134| ) has scaling dimension 1/4 and thus A c ~ g 4 / 7 . 
The chiral fixed point is thus unstable in the far IR limit, and one expects that the system will 
flow to the c = 1 fixed point of the standard three-leg spin ladder j5fj. Of course, the very 
applicability of the perturbative approach to the chiral fixed point requires that A c <C m, the 
condition which can always be satisfied for sufficiently small g. Under this condition, there exists 
an intermediate but still low-energy region A c <C E <C m where the c = 2 behaviour caused 
by frustration is dominant. The physics in this region is universal and cannot be understood 
without having recourse to the chiral fixed point. At lower energies, E <C A c , the system will 
eventually cross over to the conventional critical c = 1 behaviour. 



B. Effect of a weak interaction between surface chains at the chiral fixed point 



In this subsection, we investigate the stability of the chiral fixed point against a weak inter- 
action between the spins of the surface chains: 



01 



912S1 • S 2 gi2 (ni • n 2 + Ji ■ J 2 



(135) 



Such a term accounts for periodic transverse boundary condition in the 3-chain cylinder model 



and may change drastically the physics of the chiral fixed point. In Eq. (135), the coupling 
constant 1712 is independent of the coupling g of the current-current interaction of the model so 
that on can always consider O12 as a weak perturbation at the chiral fixed point: | <7x2 1 <C 1. 

Let us first analyse the ni • ri2 contribution at the chiral fixed point. Using the bosonizcd 
form (pq) of the staggered magnetizations of the surface chains and the Toulouse basis (100), 
one finds: 

A 2 



ni • n 2 



2n 2 a 2 



Amp^j + 2 cos (Vi^-J - cos (^/8t:^i - \/47r$; 



(136) 



Using the refermionization formula (|8J, [8^, |9^), this operator can be expressed in terms of the 
different fields occuring in the Toulouse point solution: 

iA 2 



ni • n 2 



27ran 



- (VRVL + CbCl) cos f -s 
(CrVl - tirCl) sin W< 



'87T$l 

;7r$i 



(137) 
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Since the Majorana fermions £ 3 and 77 are short-ranged fields , we can drop their contribution 
so that: 

iA 2 



ni • n 2 



27rao 



(138) 



The two terms in this equation are of different nature. One the one hand, the first contribution 
is a relevant perturbation (with scaling dimension 1) which identifies with the energy density 
operator of the Ising model corresponding to the Majorana fermion £° and drives the latter 



model out-of criticality. On the other hand, the second operator in Eq. (138) is a naively 
irrelevant contribution with scaling dimension 3 but, as we shall see below, it is crucial to 
keep it since it couples degrees of freedom of different nature: Magnetic ones described by 
the bosonic field $1 and a half of the pseudo-charge degrees of freedom stemming from the 
Majorana fermion Cr,l- 



The same analysis can be made for the current-current contribution in Eq. (135). In partic- 
ular, using the results J88| , |90| ), the leading part of this perturbation at the chiral fixed point 
affects mostly the magnetic degrees of freedom: 

-2 

Ji • J2 — 2J&JI - cos (Vi^i) (139) 



2irao 



which can be rewritten at the SU(2) point (7^ 
form: 



1, see Appendix A) in a rotationally invariant 



Ji • J 2 ~ 2Jr ■ Jjj. 



(140) 



With the two results (138, 140) at hand, we are now in position to investigate the stability 
of the chiral fixed point upon switching on a weak interaction between the surface chains. The 



effective interaction (135) at the chiral fixed point writes: 



0\2 2g 12 Jn ■ Jl + 



.9i2iA 2 
27rao 



3^£-CflCLCOs(V87T$i 



(141) 



Due to the first term, the effect of the operator (O12) on the chiral fixed point strongly depends 
on the sign of the coupling 312. 

For an antiferromagnetic interaction (i.e. g 12 > 0), we expect the succession of two transitions 
in that case. A first one (Ising transition), as soon as the coupling gi2 is switched on, the 
model is still critical but with a smaller central charge c = 3/2 since, as already emphasized, 
the Majorana fermion £° acquires a mass. As a consequence, the Z2 symmetry with respect 
to the interchange of the two surface chains is now broken. When g\2 increases, there is 
a critical value when 2512 exceeds the coupling constant 7 < of the marginal irrelevant 
current-current interaction of the effective S=l/2 Heisenberg chain associated with the bosonic 
field $1. In that case, the magnetic sector becomes massive and enters in a dimerized phase 
as in the J1-J2 problem Q. Moreover, in this massive phase, the bosonic field $1 becomes 
locked and the operator cos (-\/87r$i) acquires a non-zero expectation value: (cos (V87H&1)) ^ 
0. As a consequence, the magnetic excitation affects indirectly the pseudo-charge degrees of 
freedom by giving a mass term for the Majorana fermion £ due to the third term in Eq. 



(141): — 5 12 iA 2 (cos (v / 87r < I'i))Ci?,Ci/(2 7ra o)- Therefore, in the dimerized phase, all degrees of 
freedom are massive. In summary, the chiral fixed point is unstable against the presence 
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of an antiferromagnetic interaction between the surface chains. There are two distinct Ising 
transitions in the pseudo-charge sector whereas the spin part spontaneously dimerized for a 
sufficiently strong value of the interaction gi2. 

For a ferromagnetic interaction (512 < 0), the Majorana fermion £° still acquires a mass but 
now the bosonic field $1 and the Majorana fermion £ remain massless excitations so that the 
resulting model has central charge c = 3/2. One should note that a pseudo-charge fractionaliza- 
tion occurs in this c — 3/2 phase: A half of the degrees of freedom, described by the Majorana 
fermion £ , in the pseudo-charge sector becomes massive by an Ising transition whereas the 
other (i.e. the Majorana fermion £) is still critical. A similar phenomenon for the charge 
degrees of freedom of a one dimensional model has been reported in Ref. |}2). The fraction- 
alization manifests itself in the spin-spin correlation functions in this phase. The asymptotics 



of the correlation between spins of the central chain are still given by Eqs. (118, 119). The 



only modifications occur in the calculation of the correlation functions for the surface chains. 



For a ferromagnetic interaction (gi2 > 0), one observes from Eq. (141) that the mass Too of 



the Majorana fermion £° is positive (mo = — 3A 2 <?i2/(27rao)) so that the corresponding Ising 



model is in its disordered phase: (co) = and (fi®) 7^ 0. Using Eqs. (122, 122), one deduces 



the following estimate for the leading asymptotics of the correlation functions: 



(S* (x,t)S§ (0,0)) — 



1/1 1 



8tt 2 \(x + iuit)' 2 (x — iu±r) J 
+ (-l) l/a M — — , (142) 

(X 2 +U 2 T 2 ) 1/2 (X 2 + U 2 T 2 ) 1/8 



(#(s,r)S£-(0,0))~-^' ' ' 



(x + iu\t) (x — iu\tY 



+2(_l)*/«M — — (143) 

(x 1 + u\t 2 ) 1/2 {x 2 + u 2 t 2 ) 1/S 



where A is a non-universal constant. In this c = 3/2 phase, the slowest correlations still 
correspond to the correlator of the staggered magnetization of the surface chains and the low- 
temperature dependence of the NMR relaxation rate is now given by: 



Ti 

in contrast with the estimate (|130|) obtained at the chiral fixed point 



T (144) 



VII. CONCLUDING REMARKS 

In this work, we have studied some specific models which display in the long distance limit 
some CSL physics and belongs to a new non-Fermi-Liquid class of fixed points describing the 
IR behaviour of one-dimensional interacting chiral fermions. The possible realizations of this 
CSL state, proposed in this paper, belong to several topics of the one-dimensional Quantum 
Magnetism: the (N+l)-chain cylinder as a special frustrated spin ladder, its asymmetric doped 
version with a doped central chain as an example of a Luttinger liquid in an active environment, 
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and finally the Kondo-Heisenberg model with N channels away from half-filling as a generalized 
multichannel Kondo lattice. 

The common feature of all these models stems from the fact that their Hamiltonians in the 
continuum limit have the same structure consisting of two marginally coupled SU(2) n and 
SU(2)i WZNW models with only current-current interactions. In particular, this Hamiltonian 
decomposes into two commuting chirally asymmetric parts which transform into each other 
under the time-reversal symmetry. Each part flows to an intermediate fixed point where the 
system as a whole displays critical properties characterized by a non-trivial symmetry group 
belonging to the universality class of the CSL state |30|. All the possible lattice realizations 
are expected to exhibit critical properties in the IR limit with, in general, non-integer central 
charge. 

For the N — 2 case, by allowing anisotropic interactions in spin space, the model of two 
marginally coupled WZNW models can be solved using a Toulouse point approach. It provides 
a direct reading of the spectrum by identifying the nature of the massive and massless degrees of 
freedom in the IR limit. Moreover, the Toulouse limit solution captures all universal properties 
of the model including the SU(2) symmetric case. At the strong coupling fixed point, the 
massless degrees of freedom with central charge c = 2 consist of an effective S=l/2 Heisenberg 
spin chain and two critical Ising models which act as pseudo-charge degrees of freedom. The 
special nature of the chiral fixed point identified by the Toulouse point approach reveals itself in 
two important facts for the 3-chain cylinder model. On the one hand, the magnetic excitations 
consist of a strongly correlated bound-state made by two spinons of the surface chains and an 
antispinon of the central chain. On the other hand, the slowest spin-spin correlation functions 
corresponds to the staggered magnetizations contribution of the surface chains with an unusual 
exponent which characterizes a new universality class in spin ladders and manifest itself in 
physical quantity such as the low-temperature behaviour of the NMR relaxation rate: 1/T\ ~ 

VT. 

The Toulouse limit solution provides also a non-perturbative basis to investigate the sta- 
bility of the chiral fixed point under the presence of several operators such as backscattering 
contributions. In particular, the interchain backscattering term of the 3-chain cylinder model 
destabilizes the chiral fixed point and the system will cross over to a fixed point with central 
charge c = 1 which is expected to be the one for the non-frustrated 3-leg spin ladder with 
open transverse boundary conditions. However, there still exists an intermediate low-energy 
region governed by the chiral fixed point. This fixed point becomes also unstable upon switch- 
ing on a weak interaction between the surface chains. For a ferromagnetic interaction, the 
system exhibits an interesting critical phase with a smaller central charge c = 3/2 where a 
half of the pseudo charge degrees of freedom becomes massive by an Ising transition and the 
low temperature behaviour of the NMR relaxation rate is now: 1/Ti - T 1 / 4 . In the case of 
antiferromagnetic interaction (frustrated periodic transverse boundary conditions), we expect 
the succession of two transitions with an intermediate c = 3/2 phase followed by a massive 
phase where the magnetic sector is spontaneously dimerized. 

Regarding perspectives, it is clearly of interest to further explore the phase diagram of the 
(N+l)-chain cylinder model in the vinicity of the chiral fixed point. The experience gained 
by the Toulouse approach for N = 2 leads us to expect, though it requires a proof, that the 
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intermediate fixed point will become unstable and the system may cross over to a fixed point 
presumably characterized by a central charge c — 1. A Toulouse limit solution for the general N 
is thus clearly desirable to shed light on this problem. Using a parafermionic description of the 
SU(2)jv spin current p§( ], one can readily find a Toulouse limit for every N with the decoupling 
of a bosonic field corresponding to the bosonic held $ x in the N — 2 case. Unfortunately, the 
other critical and massive degrees of freedom strongly interact with each other: The resulting 
effective theory cannot mapped onto Majorana fermions as in N — 2 case. Otherwisely stated, 
when N exceeds two, the Toulouse limit does not correspond to a theory of free particles. The 
simplest case is the N — 4 case which is currently under study. With this respect, the situation 
is in close parallel to the multichannel Kondo problem since when the number of channels 
exceeds two, the Toulouse limit is described by a non-trivial held theory 

This analogy leads us to expect that the simple model of a Luttinger liquid in an antiferro- 
magnetic environment considered in this work and the multichannel Kondo-Hcisenberg model 
away from half-hlling should exhibit a non-Fermi-liquid low-temperature behaviour with en- 
hanced composite pairing correlations (odd-frequency pairing) as in the two channel Kondo 
problem (H^] or the one-dimensional Kondo lattice |5l],[32| . This question in the N = 2 case will 
be addressed by the Toulouse point approach in a forthcoming publication j36|. We hope that 
the chirally stabilized liquid state with all its physical properties will be observed in numerical 
simulations and in further experiments on spin ladder systems. 
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Note added 

When this work was completed, we became aware of a work by N. Andrei and E. Orignac |65[ 
who have also determined the low energy theory of the one dimensional multichannel Kondo- 
Heisenberg lattice away from half filling and found that composite order parameters induce the 
dominant instabilities for a number of channels less than 4. 

APPENDIX: A ESTIMATION OF lM AND j M 

In this Appendix, we shall compute the constants 7 m >7m that appear in the representation 



(104) of the different spin current fields in terms of the current J of the effective S=I/2 



Heisenberg spin chain at the Toulouse point. 

1. Estimation of 7 

Let us first consider the following vacuum expectation values: 

7_r = (exp (±i%/l67ri> 2 /?J ) 



7l = (exp (±iVmr$2L)). (AI) 
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In the following, we shall show that Jr = Jl = 7m, 7m being real and expresses in terms of the 
mass to of the massive modes at the Toulouse point. 

Using the definition (|59|) of the chiral bosonic field $2R,l, we have: 



T]B.,L = —F= ■ COS V167r$2i?.,L 



Cr,- 



± 



^ :sin 



(A2) 



where we have used for commodity the normal ordered form of an exponential of chiral bosonic 
fields: 



: exp (ia$ 2 R,L) 



exp (ia$ 2 fl,L) 



(A3) 



The next step is to express the previous exponential bosonic fields with scaling dimensions 2 
in terms of the different Majorana fields. This can be done by considering the following OPE 



from (A.2): 



2tt 



2tt 



vr (?) m M 



COS ( V 167T$2i?, 



1 



2n (z — lo) 



Cr (z) (r (w) 
cos ^Vl67r$2i?J : 



1 



2n (z — lo) 



4(«-w) : (a$ 2 fl) S 



(A4) 



where the fields in the second part of the previous equations are only function of lo. On the 
other hand, since rj, ( are Majorana fermions, they satisfy the following OPEs: 

f 



m (z) VR M 



- (z - co) : rjRdrjR 



2n (z — lo) 

Cr (z) Cr (w) ~ 7—7^ — — - (z-ui) ; CrQCr 
Zir (z — lo) 



(A5) 



Comparing the two short distance expansions (JAJ, A5), we end with the equivalence: 

cos (yi6TT$ 2 R) = aim (-rj R d x r} R + GAO?) , (A6) 

since for instance, in our convention, 8Cr = id x CR (z — uqj + ix). One can do the same 
computation for the left sector to find: 



cos (^/16tt$2l) = al'nr (r) L d x j] L - CzACl) 
In the same way, considering the OPE tir{z)Qr{lo) , we get: 

sin (Vl6iT$2R,L ) = ±2mal d x 7] RfL CR,L- 



(A7) 



(A8) 



This proves that 7r,l is real since at the chiral fixed point £ is a free and massless Majorana 
fermion so that (sin (V16^<£>2R,l)) = 0. Therefore, using Eqs. ( |A6| , |A7[ ) and the definition 
we obtain the following correspondence: 



1r,l = Tao i7r {(VR,Ld x i]R, L ) - (Cr,l9 x Cr,l)) 



(A9) 
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The contribution of the massless Majorana fermion £ is easy to find: 



dk f du> ike [u0+ 



=F— / dkk = =F^- 



(A10) 



where we have used a simple cut-off regularization to cure the UV singularity. To estimate 
(VR,Ld x r]R^}, we need to know the Green's function of the Majorana fermion r\ that comes 
from the massive part of the Hamiltonian at the Toulouse point describing the hybridization 
between the Majorana fermions £ 3 and r\: 



H. 



1U2 „ im / ~ \ 



(Vr -w) f] R (k, u)) (rjR {-k, -uj) £| (k, uj)) 

ie L (-k, -w) V R (k w)) (& (-k, -u) a (k, lo)) 



The Green's functions corresponding to Ti + are given by 

G+ (k,0j): 

with 

G+ 1 (k,u) = 

whereas the ones associated with Ti^ read as follows 
G_ (fc,w) 



iu; — u 2 k im 
— im iu> + v\k 



(VL i~k, -UJ) £| (fc, W)) (?/L (-fc, -W) 7/l (fc, w)) 



with 



GI 1 (k,w) 



iw — Wife — im 
+im iw + M2A; 



The Majorana Green's function for the rj field is thus given by: 

iui ± Wife 



-w) ?7fl„L (fe, W)) 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



(A16) 



(wTiki) +S 2 (fc) 

where E 2 (k) = w 2 £; 2 + m 2 and = (up + v%)/3. Integrating over the frequency, we obtain the 
following estimate: 



(VR,Ld x ij R:L ) = J dk 



E(k) 

iA 2 im 2 , / Au* 
- T-j— ± 1 o In 

47T 47TM^ \ m 

subsituting ( |A17| ) and ( |A10| ) into JA9| ), we find the value of 7 m at the Toulouse point: 

'Am* 



m -Oq 
7_r,l ~ . „ In 



,2 2 

^0 



4u 2 



(A17) 



(A18) 
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2. Estimation of 7 

The next vacuum expectation values to compute is 

1r,l = '^ao(^R,LVL,R)- (A19) 
Using the previous expression of the fermionic Green's function G±, we have: 

(&, L (-k,-u>)T, L , R (k,u)) = —. — Z (A20) 

Integrating over to, the value of 7k,£ can easily be determined: 

mao f A dk 



and thus 



214* \ m 

Let us conclude this appendix by discussing the SU(2) invariance of the spin-spin correlation 
functions between the central spins and also between the surface chains at the chiral fixed point. 
Using the Toulouse point value (|56|): gjj = Aitu* and expressing the values of 7 m ,7 m in terms 
of the ratio R = g±/g\\, we have: 

'Aoq 
2R 

l m ~Rin(^y (A23) 

In the continuum limit, the value of the physical mass is fixed. This gives the dependence of g± 
in terms of the short distance cut-off aQ. The ratio of the couplings g±/g\\ is also fixed to some 
physical value R which is a RG invariant flow of the equation (^J) . The two physical parameters 
emerging from the renormalized theory are precisely m and R apart from the different velocities 
of the problem. At the SU(2) point when k = 1, we must have 7^ = 7m = 1 so that the 



t"*-—J. eWv < A21 > 



7„ = 7 M -^lnf^|. (A22) 



R 2 In 



correlation functions of the central chain (118, 119) are rotational invariant. Using Eq. (A23) 
this fixes the value of the product Aao to 2. 



APPENDIX: B BASIC FACTS ABOUT THE BOSOZINATION OF TWO ISING 

MODELS 

In this Appendix, we briefly review some basic materials on the relation between two Ising 
models and a free bosonic field in order to fix notations for the computation of the spin-spin 
correlation function between spins of the surface chains (Section V). 

It is well known that a theory of free massless Majorana fermion with central charge c = 1/2 
describes the long distance properties of the critical two-dimensional Ising model. Two Ising 
models can be mapped onto the Gaussian model with central charge c = 1 (see for a review 
chapter 12 of Ref. |p9[ ). This procedure of doubling the Ising models has been very fruitful for 
the computation of the correlation function of the Ising model |5^- |59| , |6l|] , for the problem of 
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random impurities in the Ising model IpCf , and also for the calculation of the structure factor 
of the two-leg spin ladder JHJ] . This mapping stems from the fact that the two Majorana fields 
say 11*1,2 can be combined into a single Dirac fermion x which can be in turn be expressed in 
terms of a bosonic field $ by Abelian bosonization: 

XR V2 

with $i] = i/A to insure the anticommutation relation between xr and xl- The Majorana 
fermions ^1,2 can therefore be expressed directly in terms the $ field as follows: 



?/>l_R = —== COS f%/47T$/f 

ipiL = 7= cos (y/^n&L 
4>2R = 1 sin (VAtt^ r 



/ira 

il>2L = = sin ( V4tt$ l ). (B2) 



/7ra 



The energy density operator e a = iipaRipaL, & — 1,2, which is a mass term for the Majorana 
fermion and drives the model out of criticality, can also be written down as a function of the 
bosonic field <f> and its dual 9: 



cos ^V47r$J = i7ra (ipiRipiL + ^2Rtp2L) 

COsiV^K®) = mao (-IplRt/jlL + 1p2Rtlj2L) , (B3) 



and one can also relate all bilinears involving V'l-R.i and tp2R,L to the bosonic field $ and its 
dual 6: 

sin (y^r^j = ma (ip2RipiL ~ i>\R^2h) 

sin ^V47re^ = i7ra ("02fl^iL + ^\R^2h) , (B4) 

and 

d x <& = iy/n (lplRlp 2 R + 1plLlp2L) 

d x Q = iV7r (-ipiRipaR + ipiLifai.) • (B5) 

The Ising model labelled 1, for instance, contains apart from the Majorana fields i^ir } l that have 
conformal dimensions (1/2,0) and (0,1/2) and the energy density operator, two other fields, 
with conformal dimensions (1/8, 1/8): the order and disorder operators o\ and These fields 



are non local when expressed in terms of the Majorana fermion tJjir^l 57 5qj61|| . The Kramers- 
Wannier duality transformation maps the order operator a\ onto the disorder field In the 
ordered phase (T < T c ), (<ti) ^ 0, (^1) = whereas in the disorder phase (T > T c ), we have 
(<7i) = 0, (/ii) 7^ 0. At T = T c , both fields have a zero vacuum expectation value. At T = T c , 
the products <7i<7 2 , &1H2, Mi/^2 and \i\Oi have scaling dimension 1/4 and can be related to the 
bosonic exponentials exp(±iy / 7r<i>), exp(±i v / 7r8) f^M : 
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Hl/J,2 ~ COS (a/7T$) 

<j\02 ~ sin (\fK&) 
a\\ii ~ cos (\/7r0) 

/ii(T2 ~ sin ("v/tt©) ; (B6) 

this correspondence should also hold at small deviations from criticality. These operators are 
very useful for pratical computation since the correlation functions of the order and disorder 
fields are known exactly even out of criticality. For T = T c , we have the following power law in 
the long distance- long time limit: 

(ji(x,t) n(0,0)) = (a (x,r)a(0,0)) 

~ (3.2 +v 2 T 2y/$> ( B? ) 

v being the spin velocity of the Majorana fermion. The description of spin chains in terms of 
order-disorder parameters of the underlying Ising models has been used by Shelton et al. ||,| 
for the two-leg spin ladder and by Allen and Senechal 0| for the two-leg zigzag ladder in the 
absence of the twist term |3j|. This approach allows a simple description of the excitations of 
the system and is very useful for the computation of the correlation functions and the dynamic 
structure factors in a massive phase. 

In our problem, we shall use this approach for the computation of the spin-spin correlation 
functions of spins of the surface chains and for the study of the stability of the chiral fixed point 
perturbed by some operators. One has to express, at the Toulouse point, all the order-disorder 
operators of the underlying Ising models. First of all, the bosonic field (f- — (ipi — i/? 2 )/-\/2 
is not affected by the canonical transformation ([50], and is built from the two Majorana 
fermions £°,£ 3 . The two corresponding Ising models are labelled "0" and "3" respectively and 



using Eq. (fBfJ) with the identification (KJ, 85), we have 



jU 3 /x ~ cos (Vnifi-) 
CT 3 cr ~ sin (y/ncp-) 
cr 3 ^ ~ cos (V7rtf_) 

M 3CTo~sin(V^?-). (B8) 

A second couple of Ising models (noted "4" and "5") stems from the two Majorana fermions 
£,77 respectively and are associated with the $2 bosonic field (see Eq. (p9|)). In that case, the 



correspondence (B6) gives 

~ cos (V7T$ 2 ) 

0-50-4 ~ sin (y/TT^ 2 ) 
a 5 ^i i ~ cos (Vtt02) 

/X504 — sin (V7F02) • (B9) 

Since the £° field is decoupled from the interaction from the beginning and the Majorana 
fermion £ is massless at the Toulouse point, the Ising models noted "0" and "4" decouple from 
the other and remains critical. The correlation function of their order-disorder operators are 
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therefore given by Eq. (B7). However, the two Ising models "3" and "5" strongly interact with 
each other due to the hybridization of the Majorana £ 3 and 77 fields at the Toulouse point (see 
Eq. (64)). As shown by the spectra ( |65| , p6| ) , these two Ising models are non-critical and the 
order-disorder operators have short-ranged correlation functions characterized by the mass m 
of the massive modes at the Toulouse point. 



APPENDIX: C HYBRIDIZED ISING MODELS 

The aim of this appendix is to describe more precisely the consequence of the hybridization 
between the Majorana fermions £ 3 and r\ on their corresponding order-disorder (173,5, M3,s) Ising 
operators. In particular, we shall prove, in the following, that at the Toulouse point, one has 

(^3^5) = (0305) = M 

(ws) = (<J 3f i 5 ) = 0. (CI) 

Let us begin by recalling the effective Hamiltonian associated with the massive modes at the 
Toulouse point: 

T-Lhyb = —± {f R d x f R - £, 3 L d x £,l) - -y {md^R - rj L d x 7] L ) 

+im (CW - VR&) • (C2) 
For the computation of ground-state expectation values, we notice that the velocity anisotropy 



present in the Hamiltonian (C2) can be ignored for several reasons. On the one hand, we have 
two phcnomenological parameters in our approach: v\ and vq. One can always choose vq — v\/2 
so that all spin velocities Ui, 1*2 (pq) of the modes remain positive and v% = U2- We expect that 



the results for the vacuum expectation values (CI) obtained with this particular fine-tuning 
represent the generic situation. On the other hand, one can justify further this proposition by 
studying the structure of the equal-time fcrmionic Green's functions of the Majorana fermions. 



The nice thing is that the equal-time Green's functions corresponding to the Hamiltonian (C2) 
are in fact equal to that of a model with a single spin velocity (the average spin velocity) : 

+im (e R VL - VR&) (C3) 

where is the average of the velocities v\ and U2'- w* = {v\ + U2)/2 = (v\ + vq)/3. The 
very reason for this equivalence on the equal-time Green's functions stems from the fact that 



both Hamiltonian (C2, |C3|) have the same eigenvectors although different spectra. Therefore, 



there is a one-to-one correspondence between the static ground-state correlation functions of the 



model (C2) and that of the Hamiltonian (|C3|). The velocity anisotropy will be important when 



considering dynamical properties. Consequently, our problem is reduced to the computation 



of the ground-state expectation values (CI) with the Hamiltonian (|C3|). To this end, let us 



introduce a bosonic field $ and its dual 8: 

1 . 

„1V47t<E>j 



v/2 v^vrao 



e 



= _L= e- 1 ^^ . (C4) 
V2 v27rao 
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Using the results (B6) of the Appendix B, the different products of order and disorder fsing 
operators can be expressed in terms of the bosonic fields: 

(j, 3 fi 5 ~ cos (v^) 
03O5 ~ sin ^Vtt*) 
cr 3 // 5 ~ cos (a/^©) 

^ 3 cr 5 ~ sin (y^Qj . (C5) 



The Hamiltonian (03) can also be rewritten in the following bosonized form: 

^^((^) 2 + (^e) 2 )-£-sin(^ 

Using the substitution: 



(C6) 



4 

e -► e, 

Wft,j,6 identifies with the conventional form of a sine-Gordon model at 1 = Att: 



= + 0*©) )-J± m (vm 



(C7) 



(C8) 



In this model, the bosonic field $ is locked such as ($) =0 (m > 0). Consequently, we have 
(cos ^y 7 ^-^^ ) 7^ and (sin ^y^'I^ ) = whereas the dual field exponents exp (^tiy/nQj have 
zero vacuum expectation values. Using Eq. ( |C5| ) and the substitution (07), we thus obtain the 
above mentionned result (CI): 



(^3^5) = (0-30-5) = M^0 
(M305) = (03^5) = 0. 



(C9) 



APPENDIX: D STABILITY OF THE TOULOUSE POINT SOLUTION 

The consistency of the Toulouse approach to describe the universal properties of the chiral 
fixed point is investigated in this Appendix. 

We first begin by analysing the effect of the operators that we have neglected from the 
beginning in the Toulouse point approach. These operators are the current-current interaction of 
same chirality and the in-chain marginal irrelevant contribution that appears in the continuum 
limit of each S=l/2 Heisenberg spin chain (see Eq. (m)): 

C cc = Jok ■ (Jifl + J211) + R — * L 

Oic = JoR ■ JoL + +Jli? • JlL + J2R • J2L- (Dl) 

Let us first consider the current-current interaction of same chirality (O cc ). Using the 
bosonization of a SU(2)i spin current (jig) , we can express the operator (O cc ) in terms of 
the bosonic fields associated with each chain of the model: 
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V2n 



(d x ip R d x $ R + d x (p L d x $ L ) 



^ 2 q 2 ( C0S (y&K^R ~ V^T^R^ COS (y^'-P-R^j + R—> L 



(D2) 



Using the canonical transformation (|5C], |TJ), the current-current operator can then be written 
as a function of the different fields of the Toulouse point solution: 



1 



\/27r 



: o 2 (cos (Vl^-ii 



COS ( V87T$l COS ( Vl67T$27?. 



47T$ 



21 



+ cos WAny-Rj sin (V87r$i) sin (%/l67r$2fl + V^tHJ^i) + R Lj . (D3) 

The next step of the calculation is to extract the leading contribution of this operator in terms 
of the different critical fields at the chiral fixed point. Using Eqs. ( p9| , pij ), the non-zero 
expectation values (|7^, ^l|), and the hybridization ( gjj ) of the massive fields, the leading part 
of the current-current operator of same chirality (Df ) at the chiral fixed point is given by: 



o cc Sx^il^^ih + cos 



(D4) 



up to contributions that will give renormalization of mass and spin velocities. One can express 
the operator flD4| ) in terms of the physical spin current associated with the bosonic field 
l-i (see Eq. (§?])): 



Or. 



(JrJI 



im in 



(D5) 



At the SU(2) point, one has 7m7m = 1 (see Appendix A) and the expression (D5) can be 
recasted in a full rotational form: 



-4.7k ■ Jl- 



(D6) 



The operator (D6) coincides thus with the usual marginal irrelevant term of the continuum 
limit of the effective S=l/2 Heisenberg spin chain corresponding to the bosonic field 3>i. The 
contribution (D6) gives logarithmic corrections in the spin-spin correlation functions at the 
chiral fixed point. 

The second operator that we have neglected in our Toulouse solution is the marginally irrel- 
evant in-chain current-current interaction: 



Oi, 



JqR ' JoL + +JlR ' JlL + ^2R ■ J2L- 



(D7) 



We proceed in the same way as for the previous operator by first rewritting Oi C in terms of the 
different bosonic fields of each chain: 



o ic = i- [(d xV f - (d x df) + -L ((^$) 2 - (d x ef) + i- [{d x ^f - (d x ^f) 

— , \ „ cos (V8TTifi) — - cos (\/4Tr$) cos (V^nifi 



(D8) 



At the Toulouse point, this operator takes the following form: 
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o lc = i (W_) 2 - (a^_) 2 ) + J- ((0**i) a - (a^eo 2 ) + 1- ((a^) 2 - (d x e 
- -^-d x $id x $ 2 - -La^a^ - -i-j cos (V&F$0 cos (yi6^$ 2 ) 



2\ 
2) , 



47T u,q 



sin ^\/87r$ 1 ^ sin (Vl67r$ 2 ^ - 2 ^ 2 cos (v^^i) cos ^%/47r$ 2 ) cos (yinip 
sin ^V87r$i) sin ^V47rl> 2 ^ cos (VI^-) . (D9) 



O lc ~ ^-fl x f xifl,^ - (7 " + 2 2 2 7m) cos (V8^i) (DIO) 



27T 2 u,q 

Using Eqs. ( p9| , p4| ), and the non-zero expectation values |9l]), the leading contribution 
of the operator (0j c ) is at the chiral fixed point up to short-ranged parts and terms giving a 
renormalization of the spin velocities: 

and at the SU(2) symmetric point when 7^ = 7^ = 1 (see Appendix A), it can be reduced to 
the simple following form: 

O ic ~3J R -J L . (Dll) 

Since the original coupling constant of the in-chain current-current interaction (Oi C ) is negative, 
the effective operator ( Dll ) gives additional logarithm corrections of the spin-spin correlation 
at the chiral fixed point. 

The next check of the correctness of the Toulouse solution is to investigate the effect of small 
deviation of the coupling constant g\\ from its Toulouse-point value: g*^ = 47r(wo + i>i)/3. In 
that case, the Hamiltonian ( |6l| ) at the Toulouse point picks up an extra term: 

H CO r = -j=- (d x ip L d x <i> R + d x cp R d x $ L ) (D12) 

with 5<7|| = ej|| — gt. Using the canonical transformation (|5^, |5l| ), this operator transforms into: 

n cor = ((d x $ 1L ) 2 + (d x $ 1R ) 2 + (d x $2 L ) 2 + (d x $ 2R ) 2 ) 

{d x $ lL d x $2R + d x $ lR d x ®2 L ) ■ (D13) 

Neglecting the renormalization of the velocities wi, 2 , we end with using Eq. (|74|): 

Hear ~ 3 ' 5g 'l^ F i (d x $iLOtnR + O^irClvl) ■ (D14) 

Since the field 77 is massive, the expansion in 6g\\ does not introduce new IR singularities, 
implying that the long-distance behaviour of the correlation functions will not be modified 
except for velocities and mass renormalization. Therefore, we can conclude that the solution 
at the Toulouse point captures all universal properties of the chiral fixed point. 
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